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Abstract. We review the Airy processes; tlieir formulation and how they are conjectured 
to govern the large time, large distance spatial fluctuations of one dimensional random 
growth models. We also describe formulas which express the probabilities that they lie 
below a given curve as Fredholm determinants of certain boundary value operators, and 
the several applications of these formulas to variational problems involving Airy processes 
fT^ that arise in physical problems, as well as to their local behaviour. 
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1. Introduction 

1.1. Airy processes and the KPZ universality class. The Airy processes are a collec- 
tion of stochastic processes which are expected to govern the long time, large scale, spatial 
fluctuations of random growth models in the one dimensional Kardar-Parisi- Zhang (KPZ) 
universality class for wide classes of initial data. Although there is no precise definition 
of the KPZ class, it can be identified at the roughest level by the unusual t^^^ scale of 
fluctuations. It is expected to contain a large class of random growth processes, as well 
as randomly stirred one dimensional fluids, polymer chains directed in one dimension and 
fluctuating transversally in the other due to a random potential (with applications to do- 
main interfaces in disordered crystals), driven lattice gas models, reaction-diffusion models 
in two-dimensional random media (including biological models such as bacterial colonies), 
randomly forced Hamilton-Jacobi equations, etc. The model giving its name to the uni- 
versality class is the KPZ equation, which was introduced by Kardar, Parisi, and Zhang 

1 



AIRY PROCESSES AND VARIATIONAL PROBLEMS 



2 



[KPZ86] as a model of randomly growing interfaces, and is given by 

dth = -^{d^hy + Idlh + ^, 

where is Gaussian space-time white noise, E(^(t, a;)^(s, y)) = 6s=tSx=y 

A combination of non-rigorous methods (renormalization, mode-coupling, replicas) and 
mathematical breakthroughs on a few special models has led to very precise predictions 
of universal scaling exponents and exact statistical distributions describing the long time 
properties. These predictions have been repeatedly confirmed through Monte-Carlo sim- 
ulation as well as experiments; in particular, recent spectacular experiments on turbulent 
liquid crystals by Takeuchi and Sano [TSIO; TS12] have been able to even confirm some of 
the predicted fluctuation statistics in a physical system. 

The conjectural picture that has developed is that the universality class is divided into 
subuniversality classes which depend on the initial data class, but not on other details of 
the particular models. There are three classes of initial data which stand out because of 
their self-similarity properties: Dirac Sq, corresponding to curved, or droplet type initial 
data; 0, corresponding to growth off a flat substrate; and e^^^^ where B{x) is a two sided 
Brownian motion, corresponding to growth in equilibrium. As we will see later, each of 
these three classes correspond to concrete initial (or boundary) conditions for the discrete 
models in the KPZ class. In addition to these three basic initial data, there are three 
non-homogeneous subuniversality classes corresponding roughly to starting with one of the 
basic three on one side of the origin, and another on the other side. For one specific discrete 
model (last passage percolation or, equivalently, the totally asymmetric exclusion process) 
the asymptotic spatial fluctuations have been computed exactly for these six basic classes 
of initial data, and are given by the Airy processes: the three basic Airy processes, Airy2, 
Airyi and Airystat, and the crossover Airy processes Airy2_5.i, Airy2^BM and Airyi_>BM- 
Although these processes have been proved to arise as the limiting spatial fluctuations 
only for one model (and actually several others in the case of the Airy2 process), as a 
consequence of the universality conjecture for the KPZ class it is expected that the same 
should hold for the other models in the class. 

The purpose of this review is two- fold. First, we will explain in detail in the introduction 
the conjectural picture that we have just sketched from two different points of view: last 
passage percolation (or, more generally, directed random polymers) and the KPZ equation 
(or, more precisely, the stochastic heat equation). Along the way we will survey known 
results for these models. 

Our second purpose is to survey a collection of results for the Airy processes which 
express the probability that they lie below a given curve as Fredholm determinants of 
certain boundary value operators. These expressions have turned out to be very useful in 
obtaining some exact distributions through certain variational formulas, and in addition 
have allowed one to study some local properties of these processes. This will be the subject 
of Sections 2-4, 

Acknowledgments. JQ was supported by the Natural Science and Engineering Research 
Council of Canada. DR was supported by Fondecyt Grant 1120309 and Conicyt Basal- 
CMM. 

1.2. Directed random polymers and last passage percolation. 

1.2.1. Polymers. Consider the following model of a directed polymer in a random environ- 
ment. A polymer path is an up-right path vr = (7ro,7ri, . . . ) in (Z+)^ started at the origin, 
that is, vTo = (0,0) and nk — TTk-i £ {(Ij 0), (0, 1)} (see Figure la). On {1^+)'^ we place 
a collection of independent random weights jwij}^ .^q. The energy of a polymer path 
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segment tt of length TV is 

N 
k=l 

We define the weight of such a polymer path segment as 
(1.1) Wn{^) = e-'^-^^W = e!^^k=i<^^k 

for some fixed /3 > which is known as the inverse temperature of the model. Let ^m,n 
denote the set of up- right paths going from the origin to {M,N) G If we restrict 

our attention to such paths then we talk about a point-to-point polymer, defined through 
the following path measure on IIm,n- 

(1-2) Q'tNi^) = ZPomt(M,iV) ^^+^(^^ 

The normalizing constant 

is known as the point-to-point partition function. Similarly, if we consider all possible paths 
of length 2A'' then we talk about a point-to-line polymer, defined through the following path 
measure on Ufe=-Ar n ^N+k,N-k (that is, all paths of length 2N): 

(1-3) Q^W = ^u^^wM^)^ 

with the point-to-line partition function 

N 

Z""'^(iV)= ZP'''"'\N + k,N -k). 

k=-N 

A main quantity of interest in each case is the free energy, defined as the logarithm of 
the partition function. In the point-to-line case, another important quantity of interest is 
the position of the endpoint of the randomly chosen path, which we will denote by kjv- It 
is widely believed that these quantities should satisfy the scalings 

(1.4) log(ZP°'"*(Ar, AT)) ~ a2N + 62iV>^C2, 

(1.5) logiZ^^'^N)) ~ aiN + 6iiV>^Ci, 

(1.6) KN ~ N^T 

as N ^ oo, where the constants ai,a2 and 61,62 may depend on the distribution of the 
LOij and (3, but Ci)C2 and T should be universal up to some moment assumptions on the 
cjjj 's, and the fluctuation exponent 

X = l/3 

and wandering exponent 

e = 2/3. 

Here, and in the rest of this article, whenever we write a relation like 

as AT ^ 00, what we mean is that 

lim P ^^Vtt < ™ = IP C < m). 

V 6A^« - J _ y 



AIRY PROCESSES AND VARIATIONAL PROBLEMS 



4 



One can also have higher, d+1 dimensional versions of the model, with the paths directed 
in one dimension, and wandering in the other d. In all dimensions the scaling exponents x 
and ^ are conjectured to satisfy the KPZ scaling relation 

(1-7) X = 2C-1, 

while the universality of the limiting distributions is unclear except in d = 1. For recent 
progress on (1.7), see [Chal2; ADll; AD12]. 

Although there are few results available in the general case described above, the zero- 
temperature limit /3 — )• oo, known as last passage percolation, is very well understood, 
at least for some specific choices of the environment variables Wjj. Before introducing 
this model, we will briefly introduce the Tracy- Widom distributions from random matrix 
theory, which will, somewhat surprisingly, play an important role in the sequel. 

1.2.2. Tracy-Widom distributions. We will restrict our attention to the distributions aris- 
ing from the Gaussian Unitary Ensemble (GUE) and the Gaussian Orthogonal Ensemble 
(GOE), although these are by no means the only distributions coming from random matrix 
theory which appear in the study of models in the KPZ universality class. The reader can 
consult [Meh91; AGZIO] for good expositions on random matrix theory. 

We start with the unitary case. Let J\f{a, b) denote a Gaussian random variable with 
mean a and variance b. An N x N GUE matrix is an (complex- valued) Hermitian matrix 
A such that Aij = M{0,N/V2) + iAf{0,N/V2) for i > j and Ai^i = Af{0,N). Here we 
assume that all the Gaussian variables appearing in the different entries are independent 
(subject to the Hermitian condition). The variance normalization by N was chosen here to 
make the connection with models in the KPZ class more transparent. An alternative way 
to describe the Gaussian Unitary Ensemble is as the probability measure on the space of 
N X N Hermitian matrices A with density (with respect to the Lebesgue measure on the 
N'^ independent parameters corresponding to the real entries on the diagonal and the real 
and imaginary components of the entries above the diagonal) 

1 _J_tr42 

e 2iv"^ 

for some normalization constant Z^. If \^ , . . . , are the eigenvalues of such a matrix, 
then the Wigner semicircle law states that the empirical eigenvalue density A^~^ "l^iLi ^x'^ 
has approximately a semicircle distribution on the interval [— 2A^, 2A^]. The Tracy-Widom 
GUE distribution [TW94] arises from studying the fluctuations of the eigenvalues of a GUE 
matrix at the edge of the spectrum: if we denote by Aque(-/V) the largest eigenvalue of an 
X GUE matrix then [TW94] 

Ague(A') ~2Ar + Ari/3^2 
as — )• oo, where C2 has the GUE Tracy-Widom distribution, which is defined as follows: 

(1.8) i^GUE(s) := IP(C2 <s)= det(/ - PsKAiPs)mR), 
where i^Ai is the Airy kernel 

(1.9) KM{x,y)= / dAAi(x + A)Ai(y + A), 

JO 

Ai(-) is the Airy function. Pa denotes the projection onto the interval [a, 00), and the 
determinant means the Fredholm determinant on the Hilbert space L^(M). We will talk 
at length about the Airy kernel and related operators in later sections, so for for now 
we will postpone the discussion. Fredholm determinants can be regarded as the natural 
generalization of the usual determinant to operators on infinite dimensional spaces. We will 
review their definition and properties in Section 2. Since these determinants will appear 
often during the rest of this introduction, the reader who is not familiar with them may 
want to read Section 2 before continuing. 
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Before continuing to tlie -Fgoe we quickly note that one of the key contributions of 
Tracy and Widom [TW94] was to connect (1.8) to integrable systems. Let q{s) be the 
Hastings-McLeod solution of the Painleve II equation 

(1.10) q"{s) = 2q{sf + sq{s), 
defined by the additional boundary condition 

(1.11) q{s) ~ Ai(s) as s — ^ oo. 
Then 

Fgue(s) = e-i'r'^^(^-^)'9^W. 

The story for the Gaussian Orthogonal Ensemble is similar. An N x N GOE matrix 
is a (real-valued) symmetric matrix A such that Aij = J\f{0, N) for i > j and Ai^i = 
Af{0,^/2N), where as before we assume that all the Gaussian variables appearing in the 
different entries are independent (subject to the symmetry condition). Analogously to the 
GUE case, the Gaussian Orthogonal Ensemble can be regarded as the probability measure 
on the space ol N x N real symmetric matrices A with density 

Zn 

for some normalization constant Zf^f. As for the GUE, the Wigner semicircle law states that 
the empirical eigenvalue density for the GOE has approximately a semicircle distribution 
on the interval [— 2A^, 2A^]. The fluctuations of the spectrum at its edge now give rise to 
the Tracy- Widom GOE distribution: we denote by Aqq^(A^) the largest eigenvalue of an 
N X N GOE matrix, then [TW96] 

Agoe(^^) ~ 2iV + iVi/^Ci 
as — 7- oo, where (i has the GOE Tracy- Widom distribution, defined as 

(1.12) Fgoe(s) := P(Ci <ni) = det(/ - Po5.Po)l2(r), 
where Bs is the kernel 

(1.13) Bsix,y) = Ai{x + y + s). 

This Fredholm determinant formula for Fgoe is essentially due to [Sas05], and was proved 
in [FS05]. The original formula derived by Tracy and Widom is 

(1.14) Fgoe(s) = e-ir '^"'?WV^gue(s) 
with q as above. 

1.2.3. Last passage percolation. We come back now to our discussion about directed ran- 
dom polymers, and in particular their zero-temperature limit. We will restrict the discus- 
sion to geometric last passage percolation (LPP), where one considers a family {iOij}ijyQ 
of independent geometric random variables with parameter q (i.e. ^{uJij = k) = q{\ — q)^ 
for > 0). For convenience we also set for now cjjj- = if i or j is 0. As /3 — )• oo, the 
random path measures in (1.2) and (1.3) assign an increasingly larger mass to the path 
TT of length K > Q which maximizes the weight Wk{t^)- In the limit, the path measures 
^m'tv ^"^^ Q^N^ concentrate on the maximizing path, and the quantities which play the 
role of the free energy are the point-to-point last passage time., 

M+N 

LP°^°*(M,A^) = max V w^m 

?.=0 

and the point-to-line last passage time by 

(1.15) L'i^^(iV) = max U'°^''HN + k,N-k). 

k=-N,...,N 
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(M,iV) 




(a) (b) 

Figure 1- (a) A polymer/LPP path vr connecting the origin to 
(M, N). (b) Time s and space u axes in LPP. 



Observe that these last passage times are random, as they depend on the random environ- 
ment defined by the ujij. 

The breakthrough, which in a sense got the whole field started, was the surprising 1999 
result by Baik, Deift, and Johansson [BDJ99] which proved that the asymptotic fluctuations 
of the longest increasing subsequence of a random permutation have the Tracy- Widom GUE 
distribution. There is an intimate (and simple) connection between this model and LPP 
which we will not discuss, instead we will state the companion result by Johansson [JohOO] 
for the point-to-point LPP case: 

(L16) LP°'"*(iV, N) ~ ciN + C2A^^/^C2, 

where ci , C2 are some explicit constants which depend only on q and can be found in [JohOO] 
and ^2 has the Tracy- Widom GUE distribution. A similar result holds for point-to-line 
LPP. The longest increasing subsequence version goes back to Baik and Rains [BROl], 
while the analogue for LPP which we state here was first proved in [BFPS07] (see also 
[Sas05; BFP07]): 

(L17) L""^(A^) ~ c[N + C2N^/^Ci, 

where Ci now has the GOE Tracy- Widom distribution. 

The reason why these exact results (and others we will discuss below) can be obtained for 
geometric last passage percolation and other related models is that LPP has an extremely 
rich algebraic structure which allows one to write explicit formulas for the distribution 
of the last passage times. The algebraic structure arises from regarding the model as a 
randomly growing Young tableau, where the cell is added at time LP°™*(i,j). This 
shift of perspective relates the problem to the representation theory of the symmetric group, 
and in particular to the Robinson-Schensted-Knuth (RSK) correspondence, which is the 
main combinatorial tool used in [JohOO] to prove the following remarkable formula: 

P(LP°-'(Af, N)<s) = det(/ - P,i^^-^P,)i2(K) 

for M < N, where the Meixner kernel Kj^'^^^ is given by 

KN-l x-y 

w{x) = (*^~^'''^)) and the functions pn{x) are the normalized Meixner polynomials, i.e., 
the normalized family of discrete orthogonal polynomials pn{x) with respect to the weight 
w{x), with pn{x) of degree and leading coefficient km- A non-trivial asymptotic analysis 
of this kernel allowed Johansson to deduce that the above Predholm determinant converges 
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as — )• oo to the Fredholm determinant appearing in the definition (1.8) of the Tracy- 
Widom GUE distribution. A more detailed discussion of these facts is beyond the scope 
of this review; what the reader should keep from this discussion is that the exact results 
which we are discussing depend crucially on what is usually referred to as exact solvability 
or integrability: the availabilty of (extremely non-trivial) exact formulas for quantities of 
interest. These formulas arise from the very rich algebraic structure present in some (but 
by no means all) models in the KPZ class. For a recent survey on this subject see [BG12]. 

As part of the general KPZ universality conjecture, one expects that (1.16) and (1.17) 
hold not only for LPP, but in general for any /3 > 0. In other words, the belief is that 
in (1.4) and (1.5), the random variables C2 and have respectively the Tracy-Widom 
GUE and GOE distributions. There has been only partial progress in proving this con- 
jecture for point-to-point directed polymers (and virtually none in the point-to-line case), 
the difficulty lying in the lack of exact solvability. Versions of this conjecture have been 
proved for two related models in the point-to-point case: the continuum random polymer 
in [ACQll] (building on results of [TW08a; TWOSb; TW09]) and the semi-discrete poly- 
mer of O'Connell and Yor in [BCll; BCF12] (see also [0'C12]). In the setting of discrete 
directed random polymers, [COSZll] showed that if the weights are chosen so that —Wij is 
distributed as the logarithm of a Gamma random variable with parameter 9i + 6j (for some 
fixed 9iS and ^j's) then the model is exactly solvable in the sense explained above. This 
was later used in [BCR12b] to prove that the asymptotic fluctuations of the free energy 
of the point-to-point polymer (at least for low enough temperature) have the conjectured 
Tracy-Widom GUE distribution. 

1.2.4. Spatial fluctuations and the Airy processes. The Airy processes arise from LPP when 
we look not only at the fluctuations of the free energy at a single site, but instead at several 
sites. To this end, we deflne the rescaled point-to-point process u 1— )• H^^^^{u) by linearly 
interpolating the values given by scaling LP°'°*(iVf, A^) through the relation 

(1.18) LP°"'\N + u,N -u) = ciN + C2N^/^Hf;'''\c3N~^/\) 

for u = —N, . . . ,N, where the constants q have explicit expressions which depend only 
on q and can be found in [Joli03]. Observe that this corresponds to looking at the free 
energy at a line of slope —1 passing through {N,N). The limiting behavior of is 
described by the Airy2 process A2- This process was introduced by Prahofer and Spohn 
[PS02], and is deflned through its flnite-dimensional distributions, which are given by a 
Fredholm determinant formula: given xq, ■ ■ ■ ,Xn S and ui < • • • < in M, 

(1.19) nMui) <Xi,.. .,A2{Un) < Xn) = det(/ - f^/'i^lf f^/2)L2(|„^,...,„„|xR), 

where we have counting measure on {ui, . . . , Un} and Lebesgue measure on M, f is deflned 
on {m, . . . ,Un} X M by 

(1-20) f('"j,a;) = lxe(xj,oo), 

and the extended Airy kernel [PS02; FNH99; Mac94] is defined by 

(1.21) K^.An,^,u,^}-<^_ ^.^^ ^ ^ ^ ^ 

Although it is not obvious from the definition, the Airy2 process is stationary (this will 
become clear in Section 3.1), and as should be expected from (1.16), F{A2{u) < m) = 
FQUE{m) for all u. There is a close connection, which we will explain in Section 1.5, 
between the Airy kernel K^i appearing in the definition (1.8) of the Tracy-Widom GUE 
distribution and the extended kernel K'^^. 

The precise result linking the point-to-point LPP spatial fluctuations to the Airy2 process 
is due to Johansson (see also [PS02]): 
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Theorem 1.1 (Johansson [Joh03]). There is a continuous version of A2, and 

//P°-t(^) yA2{u)-u^ 

in distribution in the topology of uniform convergence of continuous functions on compact 
sets. 

In the LPP picture, the "time" variable (which we wiU denote by s) flows in the (1, 1) 
direction of the plane, while "space" (which, as above, we will denote by u) corresponds 
to the direction (1,-1) (see Figure lb). In this sense. Theorem 1.1 describes the spatial 
fluctuations of the point-to-point last passage times as time s — )• 00. 

One can think of extending the LPP model to paths starting at s = with any space 
coordinate, i.e., paths which start at any point of the form {k, —k), k G Z. To recover 
point-to-point LPP one simply sets uJij = whenever i < or j < 0, which is easily seen 
to be equivalent (from the point of view of last passage times) to forcing our paths to start 
at the origin. In this sense, point-to-point LPP and the Airy2 process correspond to the 
5o (also known as delta, narrow wedge or curved) initial data (see Figure 2a). Note that 
in this case we only assign positive weights to sites such that s > \u\. To recover the flat, 
stationary and mixed initial data which we introduced earlier, we need to assign weights 
to sites such that s < \u\. 

Remark 1.2. The results for the flat, stationary and mixed initial data have been proved 
in settings which differ slightly from the one introduced here. To avoid additional notation 
and complications, we will state the results on the present setting, and refer the reader to 
the corresponding references for more details on the differences. 

We start with the flat initial data. It corresponds to extending the weights Uij to be 
independent geometric random variables with parameter q whenever i + j > and setting 
ujij = otherwise (see Figure 2b). This corresponds to letting our paths start at any site 
in the line s = but not attach any additional weights along that line, which explains the 
name flat. The corresponding point-to-line rescaled process may be defined as follows: first 
we extend the definition of last passage times to accomodate the flat initial data, 

2i+AI+N 

LP7*(M, TV) = max max ^ 

i££ 'reii{i,-i)^{i\/,jV) 

with self-explanatory notation, and then we define the rescaled process u 1— )• H^^'^{u) by 
linearly interpolating the values given the relation 

LlTi^ + u,N-u) = ciN + C2N^^''Hf}'%C3N-^/\) 

for u = —N,...,N. The flat initial data gives rise to the Airyi process Ai, which was 
introduced by Sasamoto [Sas05], and is defined through its finite-dimensional distributions, 

(1.22) P(^i(ni) < 6, ■ • .,Aiiun) < Cn) = det(/ - iK^'^) LH{u,,...,u„}xRr 
with f as in (1.20) and 

(1.23) Kr\u,^;u',a = -^=i^=exp(-^ff^) 1„,>„ 

^y4^^{u' -u) \ 4(ii' - u) J 

+ Ai(e + i' + {u' - uf) exp(('u' - n)(e + O + | in' — n) ) . 

The Airyi process is stationary, and as should be expected from (1.17), its marginals are 
given by the Tracy- Widom GOE distribution: ¥'[Ai{u) < m) = Fgoe(2?7t-) for all u. 

Theorem 1.3 ([BFPS07; BFP07; BFSOSa]). 

H^^'^u) > 2^3^^(2-2/3^) 

N^oo 
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(a) Narrow wedge (b) Flat 
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(c) Stationary 



Figure 2. Schematic representation of the LPP models with asymptotic 
spatial fluctuations given by: (a) Airy2; (b) Airyi; (c) Airy stat- 
in the sense of convergence of finite- dimensional distributions (on a slightly different setting 
than the one presented here, see Remark 1.2). 

The powers of 2^/^ in the above limit should be regarded as an arbitrary normalization 
(in fact, one could have defined the Airyi process as this scaled version of it). The ap- 
pearence of these factors has to do with the fact that the natural scaling in the definition 
of these quantities differs between random matrix models and models such as LPP or di- 
rected polymers (as, for instance, in (1.25) below). See Section 2 of [CQR12] for a related 
discussion. 

The stationary initial data is slightly more cumbersome to introduce. The name station- 
ary comes from the fact that for the closely related totally asymmetric exclusion process 
(TASEP), this initial condition corresponds to starting with particles placed according 
to a product Bernoulli measure with parameter 1/2, which is stationary for the process. 
Translated to LPP, this initial condition corresponds to the following. Let (5'n)nez be 
the path of a double-sided simple random walk on Z with 5*0 = 0, which we assume 
to be independent of the weights Wjj. We rotate this random walk path by an angle 
of — 7r/4 and then put it along the s = line by defining the (random) discrete curve 
7o = {(|(S'(i) + i), \ {S{i) — i)) , i G Z}. We then extend the weights Wjj to be independent 
geometric random variables with parameter q whenever (i,j) hes above 70 and = 
otherwise (see Figure 2c). The corresponding stationary rescaled process Hff''^{u) can be 
defined analogously to the previous cases, by maximizing over paths starting at 70 and 
going to the anti-diagonal line passing through {N, N). It gives rise to the Airystat process 
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-^stat • Its definition is also given in terms of finite-dimensional distributions involving Fred- 
holm determinants, but the formulas are a lot more cumbersome. We will not need the 
exact formulas, so we refer the reader to [BFPIO] for the details. Despite its name, .Astat is 
not stationary as a process. In fact, due to the connection with stationary TASEP, .Astat is 
just a standard double-sided Brownian motion, but with a non-trivial random height shift 
at the origin given by the Baik-Rains distribution, see [BROO]. The convergence result in 
this case is the following: 

Theorem 1.4 ([BFPIO]). 

HT\u) >AstM 

Af— s>oo 

in the sense of convergence of finite- dimensional distributions (on a slightly different setting 
than the one presented here, see Remark 1.2). 

The mixed initial conditions can be obtained by placing one condition on each half of 
the line « = 0. We will explain how this is done in the case of the half- flat, or wedge-^flat 
initial data, and leave the examples leading to ^2-s>bm and Ai^bm to the interested reader 
(see [BFS09; CFPIO]). To obtain the Airy2_>i process we extend the weights ujij to be 
independent geometric random variables with parameter q whenever i, j > 0, or i -|- j > 
with i < 0, setting coij = for all other sites. The half-flat rescaled process H^^^~^^^^{u) 
is obtained as in the previous cases, and gives rise to the Airy2_s.i process ^2^1- It was 
introduced by Borodin, Ferrari, and Sasamoto [BFSOSb], and is given by 



P(^2->l(ni) < 6, • • .,A2^l{Um) < U) = det{l - fi^2'^if)^2(|,^, 

with f as in (1.20) and 



.,Mm}xI 



u'>u 



If f ^w^ /3+u'w'^-(,'w 

+ / dw / dz 



(27ri)2 J^_ ^z'^i/s+uz'^-^z (^z — w){z -\- w)^ 

where = ^ — u'^lu<o, = C ~ (^^')^l«'<o and the paths 7+, 7- satisfy —7+ C 7_ with 
7+ : e''^+oo e''"^+oo, 7_ : e~'"f'-oo e'^-00 for some 0+ G (7r/3,7r/2), G (vr/2,7r - 
(p+). As could be expected from the above description, the Airy2-!.i process crosses over 
between the Airy2 and the Airyi processes in the sense that A2^i{u + v) converges to 
2^/^Ai{2~'^/^u) as v — )• 00 and A2{u) when v — )• —00. The convergence result is the 
following: 



Theorem 1.5 ([BFSOSb]). 



^half-linc(^) - uHu<0 > ^2^l(^) 

A"— !-oo 



in the sense of convergence of finite- dimensional distributions (on a slightly different setting 
than the one presented here, see Remark 1.2). 

Some of these results have been extended to the case where the points at which one 
computes the corresponding finite-dimensional distributions do not all lie in the same anti- 
diagonal line, but instead fall in certain space-like curves lying close enough to such a line, 
see e.g. [CFPIO] for more details. 

From the definitions it is clear that the basic three Airy processes A2, Ai, and .4stat are 
invariant under A{u) 1— )• A{—u), but the mixed cases are not. 

Since all initial data are superpositions of Dirac masses, there is a sense in which the 
Airy2 process is the most basic. For example, using the fact that point-to-line last passage 
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times are computed simply as the maximum of point-to-point last passage times, Johansson 
[Joh03] obtained the following celebrated formula as a corollary of (1-16) and Theorem 1.1: 

(1.25) P( sup {A2{x) - x^} < m ) = FGOE(4^/^m). 

A direct proof of this formula was later provided in [CQR12]. The argument used in this 
second proof starts with a different expression for the finite-dimensional distributions of 
A2 in terms of the Fredholm determinant of a certain boundary value operator. This type 
of formula, and their extensions to continuum statistics, are the starting point of most of 
the results we will survey in Sections 3.2 and 4. For example, as described in Section 4.2, 
they allow to compute the asymptotic distribution of kjv, the position of the endpoint in 
the maximizing path in point-to-line LPP. 

Extrapolating from (1.25) leads to a conjecture that the one-point marginals of the other 
Airy processes should be obtained through certain variational problems involving the Airy2 
process. To state the precise conjectures we turn to the stochastic heat equation, whose 
logarithm is the solution of the KPZ equation. The advantage of this model over LPP and 
other discrete models is that it is linear in the initial data, and hence the heuristics are 
more easily stated in that context. The disadvantage is that most of the argument relies 
on conjectures based on universality. 

1.3. The continuum random polymer and the stochastic heat equation. We now 

consider the continuum version of the finite temperature discrete random polymers (1.2) 
and (1.3). The (point-to-point) continuum random polymer is a random probability mea- 
sure Pj^'^ on continuous functions x{t) on [0,T] with x(0) = and x{T) = x with formal 
weight 

given to the path x(-), where ^(t,x), t > 0, x G M is space-time white noise, i.e. the 
distribution- valued Gaussian variable such that for smooth functions (/? of compact support 
in M+ X M, := f^^^-^dt dx ip{t, x)^{t, x) are mean zero Gaussian random variables 

with covariance structure E[{(pi,S^){ipi,^)] = {ipi, 992)- One can also think of the continuum 
random polymer as having a density 

g-/3/o^d*5(t,^W) 

with respect to the Brownian bridge. Neither prescription makes mathematical sense, 
but the second one does if one smooths out the white noise ^(t, x) in space. Removing 
the smoothing, one find that there is indeed a limiting measure supported on continuous 
functions C[0,T] which we call -P^'f - In fact, it is a Markov process, and one can define it 
directly as follows. Let z{s, x, t, y) denote the solution of the stochastic heat equation after 
time s > starting with a delta function at x, 

(1.26) dtz = \dlz - lie^z, t>s,yeR, z{s,x, s,y) = 6x{y). 

It is important that they are all using the same noise ^. Note that the stochastic heat 
equation is well-posed [Wal86] . The solutions look locally like exponential Brownian motion 
in space. They are Holder | — ^ for any S > in x and \ — ^ for any 6 > in t. In 
fact, exponential Brownian motion e^^^^ is invariant up to multiplicative constants, i.e. if 
one starts (1.26) with e^^^^ where B(x) is a two-sided Brownian motion, then there is a 
(random) C{t) so that C{t)z{t, x) is an exponential of another two-sided Brownian motion 
[BG97]. Pj-'^ is then defined to be the probability measure on continuous functions x{t) on 
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[0, r] with x(0) = and x(T) = x and finite dimensional distributions 

_ z{0, 0, ti,Xi)z{ti,Xi,t2, X2)--- z{tn-l,Xn-l,tn, Xn)z{tn, Xn,T, x) 

z{0,0,T,x) ^" 

for < ti < t2 < ■ ■ ■ < tn < T . One can check these are a.s. a consistent family of 
finite dimensional distributions. This holds basically because of the Chapman-Kolmogorov 
equation 

du z{s, X, T, u)z{t, n, t, y) = z{s, x, t, y) 

for s < T < t, which is a consequence of the linearity of the stochastic heat equation. 

Note that the construction is for each T > fixed. Unlike the usual case of diffusions, the 
measures are very inconsistent for varying T. One should imagine that the polymer paths 
are peeking into the future to see the best route, so the measure depends considerably on all 
the noise in the time interval [0, T]. We can also define the joint measure = Pj''^<^Q{C) 
where Q is the distribution of the ^, i.e. the probability measure of the white noise. 

Theorem 1.6 ([AKQ12b]). 

(i) The measures -P^'| andP^^. are well-defined (the former, Q-almost surely). 

is a Markov process supported on Holder continuous functions of exponent 
^ — (5 for any 6 > 0, for Q-almost every ^. 
(Hi) Let t"^ = Then with P^^ probability one, we have that for all < t < 1 

[2"tJ 

{x{tl)-x{tl_,) f ^t, 

fc=l 

i.e. the quadratic variation exists, and coincides with the one obtained for Fj^ ^ (the 
Brownian bridge measure), 
(iv) Pj''^ is singular with respect to P^^ (the Brownian bridge measure) for Q-almost 
every ^. 

So the continuum random polymer looks locally like, but is singular with respect to, 
Brownian motion. One can also define the point-to- line continuum random polymer Pj,, in 
the same way as in the discrete case. For large T, one expects Varp/3 (x(T)) ~ T^/^ in the 

point-to-line case or Varp/3 {x{T/2)) ~ T'^/^ in the point-to-point case. Here the variance 

is over the random background as well as P'^'^- The conditional variance given ^ should be 
much smaller. 

If z{t,x) is the solution of (1.26) then h(t,x) = — (3~^ log z{t,x) can be thought of as 
either the (renormalized) free energy of the point-to-point continuum random polymer, or 
the Hopf-Cole solution of the Kardar-Parisi- Zhang equation, 

(1.27) dth = -^{d^hf + ldlh + ^, 

for random interface growth. Since log z{t,x) looks locally like Brownian motion, (1.27) 
is not well-posed (see [Haill] for recent progress on this question.) If ^ were smooth, 
then the Hopf-Cole transformation takes (1.26) to (1.27). For white noise ^, we take 
h{t, x) = log z{t, x) with z{t, x) a solution of (1.26) to be the definition of the solution 

of (1.27). It is known [BG97] that these are the solutions one obtains if one smooths 
the noise, solves the equation, and takes a limit as the smoothing is removed (and after 
subtraction of a diverging constant). They are also the solutions obtained as the limit 
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of discrete models like asymmetric exclusion in the weakly asymmetric limit [BG97], or 
directed polymers in the intermediate disorder limit [AKQ12a]. 

To understand the intermediate disorder limit we consider how the KPZ equation (1.27) 
rescales. Let 

Kit,x) = e'^hie-H^e-^x) 
Recall the white noise has the distributional scale invariance 

i{t,x)^^'e'^i{e%e^x). 

Hence, setting (3 = 1 for clarity, 

Because the paths of h are locally Brownian in x we are forced to take a = 1/2 to see 
non-trivial limiting behaviour. This forces us to take 

z = 3/2 

The non-trivial limiting behaviour of models in the KPZ universality class are all obtained 
in this scale. 

On the other hand, if we started with KPZ with noise of order e^/^, 

dth = -\{d^hf + \dlh + e^/^i, 

then a diffusive scaling, 

h^{t,x) = h{e^^t, e^^x), 

would bring us back to the standard KPZ equation (1.27). This is the intermediate disorder 
scaling in which KPZ and the continuum random polymer can be obtained from discrete 
directed polymers. It tells us that if we set 

in (1.1) then the distribution of the rescaled polymer path 

Xe{t) := ex|^e-2fj < t < T 

will converge to the continuum random polymer, with temperature c/3 (see [AKQ12a] for 
details). 

1.4. General conjectural picture for the SHE. Define At from the solution of (1.26) 
by 

(1.28) z{0,y;t,x) = _i^e-^-^+2-/3,i/3^,(2-i/3i-2/3(._,))_ 

At{-) is called the crossover Airy process, the key conjecture being 

(1.29) At{x) ^ A2{x) 

This is known in the sense of one-dimensional distributions (see [ACQll], where (1.29) 
is Conjecture 1.5). A non-rigorous derivation based on a factorization approximation for 
the Bethe eigenfunctions of the (5-Bose gas can be found in [PS 11]. Note however that the 
factorization assumption is almost certainly false. 

Now one tries to use the linearity of the stochastic heat equation to solve for general 
initial data z{0,x) = zq{x), 

(1.30) z{t,x) = I ^e-^-^+2-V3,i/3^,(2-i/3,-2/3(,_,))^^^^^^ 

It is not hard to see that the equality is correct in the sense of one-dimensional distributions, 
but not more. If one wants, for example, joint distributions of z{t, Xi) for more than one Xi, 
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then one needs to enhance the crossover Airy process in (1.28) to a two parameter process 
At(2-^/3t~2/3^^2-i/3t-2/3y), xhe conjectural limit of this is a two parameter process we 
call the Airy sheet. However, we do not even have a full conjecture for its finite dimensional 
distributions, though some properties can be described (see [CQll]). 

Calling X = and y = l^^^'^t^'^^^y and starting with initial data zo{x) = 

exp{2-V3ti/3j(2-i/3i-2/3^)|^ 

we can rewrite the exponent in (1.30) as 

2-1/3^1/3 ^Atix -y)-{i- yf - /(y)] - 

so that for large t the fluctuation field 2-'^/'^*^-'^/^ [logz(t, x) + + log(\/27rt)] is well 
approximated by 

sup {A2{x -y)-{x- yf - f{y)]. 

The type of initial data would appear to be quite restrictive, but actually this picks out 
the appropriate self-similar classes. The easiest example is the flat case / = 0. We obtain 
the statement 

(1.31) Ai{x) sup{^2(y-x)-(y-x)2} 

y 

in the sense of one-dimensional distributions. Since the left hand side is just the GOE 
Tracy- Widom law this is the well known theorem of Johansson (1.25) once again. 

If one starts with a two sided Brownian motion, then the required self-similarity of this 
initial data is just the Brownian scaling and one arrives at 

Atat(a:) sw-p {A2{y - x) - {y - xf - B{y)] . 

The mixed cases require a tiny bit more care. Let's explain the heuristics first for the case 
of the Airy2_!.i process. Starting from the step initial data z{0, x) = lx>o the prediction is 

(1.32) - log z{t, x) « ^x2l^<o + ^t + log(\/2^) - 2-^/h^/^A2M2-^/h-^/^x). 
On the other hand, by linearity we have for each fixed x, in distribution, 

(1.33) z{t,x)= dyz{0,y;t,x)= ^e-^-^£+^-^'^*^'^^*(2-/3.-/3(.-,))_ 

Jo Jo 
Comparing with (1.32) we deduce that the processes supj^>o (^2(3^ — v) — {x — y)^) and 
A2~^i{x) —x'^lx<:o should have the same one-dimensional distribution or, equivalently, that 

(1.34) A2^i{x) - xHx<o = sup{^2(y) -y^} 

y<x 

for each fixed x G M. This distributional identity has actually been proved rigourously, and 
its proof is based on the methods we will survey in Sections 3 and 4 (see Theorem 4.6). 

The same heuristic argument works for the other two crossover cases. If we let z{0, x) = 
e^^^^lx>o, where B{x) is a standard Brownian motion, then (1.32) and (1.33) are replaced 
respectively by 

- log Zit, X) « ^^xHx<0 + Mt + log(\/2^) - 2-1/3^1/3^2^3^^(2-1/3^-2/3^) 

and 

z{t,x) = r dyz{0,y-t,x) = T 1 e-T-M+^(^)+2- W/3a.(2-V3,-./3(._,))^ 
Jo Jo ^ 

and now the same scaling argument allows to conjecture that 

A2^bm{x) - x^lx<o = sup{A2{y) + B{x-y)-y'^) 

y<x 
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for each fixed a; G M, where now B{y) is a Brownian motion with diffusion coefficient 2. 
An analogous argument with z{0,x) = lx<o + e^^^^lx>o translates into conjecturing that 

^i^BM(a;) = sup {A2{y) +B{x - y)ly<x - y"^) 

for each fixed x G M. As we explained, these equalities in distribution will only hold in the 
sense of one-dimensional distributions, i.e. for each fixed x. 

The strategy used in the proof of (1.34) is considerably more difficult to implement for 
the other two crossover cases (see Section 4.4 and the discussion at the end of Section 1.2 
in [QR12b]), and in fact these identities remain conjectures for now. In work in progress 
[CLW] obtain an improved version of the slow decorrelation result proved in [CFP12], which 
should allow to prove a general version of formulas for last passage times in last passage 
percolation in terms of variational problems for the Airy2 process. In particular, such a 
result would give a proof of these conjectural formulas. 

1.5. Determinantal formulas and extended kernels. As we already mentioned, the 
results we will survey in Sections 3 and 4 are based on alternative Fredholm determinant 
formulas for the finite-dimensional distributions of the Airy processes. We will introduce 
these formulas in Section 3, but before doing that let us explain why the original extended 
kernel formulas are natural. We will first do this in a simpler setting, namely random point 
processes on finite sets. Then we will explain how similar arguments can be used to derive 
the formula (1-19) for the finite-dimensional distributions of the Airy2 process. 

1.5.1. Extended kernels and the Eynard-Mehta Theorem. Let X he a. finite set. A random 
point process on is a probability measure on the family 2'^ of subsets of X, which we 
think of as point configurations. A random point process is called determinantal if there 
exists a \X\ x \X\ matrix K with rows and columns indexed by the elements of X such that 

p{A) := P({X G 2^ : A C X}) = det(if U), 

where P is the probability measure underlying the point process and K\a is the submatrix 
of K indexed by A, 

Ka=[K{x^x'%^^,^^. 

The function p is called the correlation function of the process, and K is called its corre- 
lation kernel. 

We are interested in a particular type of random point processes. Let X^, . . . , Af" be 
n disjoint finite sets. We consider a point process supported on /cn-point configurations 
with the property that there are exactly k points in each X^ . The probability of such 
a configuration is given as follows: given collections of points {2;]}j=i,...,fe ^ X^ for j = 
1, . . . , n, we set 

(1.35) P({{xi}i=i,...,fc U • • • U 

= Z-Met [Ux])\l^^,dei [wM,x])\l^^,---dei [Wn-i{xr\x^)\l^^,dei [V'^(x,")] ',,=1- 

where the (f)iS are some functions on X^, the ijjiS are some functions on X^ and the WiS 
are matrices with rows indexed by X'^ and columns indexed by X^~^^. The normalization 
constant Z is chosen so that the total mass of the measure is 1. We are assuming implicitly 
that the right hand side above is non-negative for any admissible point configuration. 

Write ^> for the k x \X^\ matrix with k rows and columns indexed by elements of X^ 
which is defined by ^i^x = (pii^) for 1 < i < k and x G X^. Similarly, write ^ for the 
X k matrix with k columns and rows indexed by X'^ which is defined by "i>x,i = i^iix) 
for 1 < i < k and x G X"'. Furthermore, define the k x k matrix 
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We will assume that det(M) 7^ 0. Under this assumption it can be shown (see e.g. [BR05]) 
that the normalization constant Z in (1.35) equals det(M). 

The Eynard-Mehta Theorem states that a random point processes defined as in (1.35) is 
determinantal. Moreover, the theorem gives an explicit formula for the correlation kernel. 
The precise statement is the following: 

Theorem 1.7 ([EM98]). The random point processes defined by (1.35) is determinantal. 
Its correlation kernel is the block matrix K with nx n blocks, such that the {i,j)-block has 
rows indexed by and columns indexed by , and is given by 

Kij = Wi--- Wn-i'ifM-^^Wi ■ ■ ■ Wj^i -Wi--- Wj-i. 

For a simple proof of this result see [BR05]. Remarkably, the inverse can be com- 
puted, or at least approximated, in many cases of interest. 

The connection with the models we have discussed so far is through certain families of 
non-intersecting paths. The Airy2 process can be obtained directly as a limit of the top line 
of several different families of non-intesecting paths, one of which is presented in Section 
1.5.2 (for some others see [Joh06]). For the other Airy processes presented in Section 1.2.4 
the connection with non-intersecting paths is less immediate (see for instance the discussion 
preceding Lemma 3.4 in [BFPS07]), but in any case enough of the above structure remains, 
and the proofs still rely crucially on a version of Theorem 1.7. On the other hand, we may 
think of the fen-point configurations where the measure defined in (1.35) is supported as 
defining a family of k (in principle not necessarily non-intersecting) paths. For example, 
the first path would be expressed by {x\,x1, . . . ,Xi). It turns out, as we will see below, 
that probability measures on families of k non-intersecting paths on X^ U • • • U X"' are 
naturally given by expressions like (1.35), and hence have a determinantal structure. If the 
sets X^ are endowed with some total order and we assume that our non-intersecting paths 
are arranged so that {x\,x1, . . . , x^) is the top path, then one can prove (see e.g. [Joh03]) 
that 

(1.36) F{x\ <zi,...,x'l< Zn) = det{I + PKP), 

where Zi G X^, K is the correlation kernel given by the Eynard-Mehta Theorem and P is 
block-diagonal matrix with n diagonal blocks defined so that, for i = 1, . . . ,n, Pi^i has rows 
and columns indexed by X'^ and is given by {Pi^iv)j = l^t^^.Vj. This should be compared 
with an expression like (1.19). 

If we go back to thinking about these paths as defining a random point process, then 
they are given by a measure on fen-point configurations on X^ U • • • U X^ . Therefore, if the 
process is determinantal, its correlation kernel necessarily has to be a matrix with rows and 
colums indexed by X^ U • • • U X^ . The Eynard-Mehta Theorem implies moreover that the 
correlation matrix is partitioned naturally into nx n blocks, with the {i,j) block having 
rows indexed by X"^ and colums indexed by X^ . To see how this structure relates with the 
extended kernels introduced in Section 1.2.4 for the Airy processes, we make the following 
observation. An operator T acting on L^{{ti, . . . , tn} x M) can be regarded as an operator- 
valued matrix (Tjj) . ^ with entries Tjj (acting on L^(M)), which acts on / G L^(M)" 

as {Tf)i = X]j=i '^i,jfj (more precisely, we are using the fact that L^({ii, . . . , t„} x M) and 
®ie{ti in}-^^^^) ^'■^ isomorphic as Hilbert spaces). Hence an extended kernel formula 
like (1.19) can be thought of as the determinant of an n x n matrix whose entries are 
operators on L^(M). Similarly, we may think of (1.36) as the determinant of an n x n 
matrix whose (i, j) entry maps L^(Af*) to L'^{X^). Since the Airy processes live on the real 
line instead of finite sets, this latter spaces are replaced by L^(M). 

In the next section we will explain how these ideas can be used to derive the extended 
kernel formula for Airy2. 
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1.5.2. Derivation of the Airy2 process from Dyson Brownian motion. The original deriva- 
tion of the Airy2 process was done in [PS02] using quantum statistical mechanical argu- 
ments, while Johansson's proof of Theorem 1.1 relies crucially on the connection between 
LPP and the Robinson-Schensted-Knuth algorithm, which provides a family of discrete 
non-intersecting paths, the top line of which converges to A2- We will briefly explain the 
derivation of the Fredholm determinant formula for the Airy2 process using a different 
model, the Dyson Brownian motion. We refer the reader to [TW07; TW04; Joh03; Joh06] 
for more details on the derivation of Airy2 from non-intesecting paths. 

Consider the evolving eigenvalues of an N x N GUE matrix with each (algebraically 
independent) entry diffusing according to a stationary Ornstein-Uhlenbeck process. We 
write the eigenvalues at time t as A^(t) = {X^ (t), . . . , X^{t)) so that Xi{t) decreases in 
i. This eigenvalue diffusion, called the stationary GUE Dyson Brownian motion, can 
be written as the solution of a certain A^-dimensional SDE, and it can be shown that 
it is stationary, with distribution given by the eigenvalue distribution of an x GUE 
matrix. Moreover, the paths followed by the A^ eigenvalues almost surely form an ensemble 
of non-intersecting curves. 

Suppose we look at this eigenvalue diffusion at times ti < ■ ■ ■ < tn, and we condition 
the A^ paths to be pairwise non-intersecting. To investigate the transitions between tm and 
tm+i, suppose we condition this eigenvalue diffusion to start at time tm at Xf (tm) = Xi for 
some fixed xi < ■ • • < xn, and we also fix destination points yi < • • • < vn- Let pt{x, y) be 
the transition probability density of a (one-dimensional) Ornstein-Uhlenbeck process from 
X at time to y at time t. Then, in this setting, the Karlin-McGregor Theorem [KM59] 
implies that the transition probability density for these A^ non-intersecting paths to end 
at the prescribed destination points yi, . . . ,y]\f is given by a constant times 



(1-37) det[pt,^^,-t,Axi,yj)] 



N 



The transition function Pt^+x-tm corresponds then to the matrix Wm in (1.35). Of course, 
since our paths take values in M now, we no longer have a matrix, but the Eynard-Mehta 
Theorem still holds in this setting (see e.g. [TW04]). The functions and tpi in (1.35) 
are related to the (stationary) marginals for A^, and in this case are equal (due to sta- 
tionarity) and expressed simply in terms of Hermite polynomials. The result after further 
computations and using (1.36) is the following [TW04]: given xi, . . . ,Xn G M, 

P(Af (tl) < XI, . . . , Af (t„) < Xn) = det(/ - fi^frm,ivf)i2„,^,...,,„},K)' 

where f is defined as in (1.20) and K^^_^^ is the extended Hermite kernel 

( N-l 

^ e^(^-*)(/p,.(x)v;fc(y) iis>t, 
-^Hrm,Af ('5) 3;; t, y) — ' ^ ^ 

■^e''^'-^'^^k{x)ipk{y) ifs<t, 

k=N 

and where ^pk{x) = e~^^^'^pk{x) with p^ the fe-th normalized Hermite polynomial (so that 
llv'felb = !)• Note that Xi is the top line of our family of non- intersecting paths, so this 
probability is the same as the probability that all paths stay below the Xj's. Note also the 
similarity between this formula and the formula (1.19) for A2- 

The kernel A'g^^ ^ has a nice algebraic structure. Writing 

D = -l{A-x'^ + l), 

(i.e. Df{x) = — |(/"(x) — (x^ — l)/(x))), the harmonic oscillator functions (pk satisfy 
Dipk = kifk, and moreover {(/5^,}^,>o forms a complete orthonormal basis of L^(M). Define 
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the Hermite kernel as 

N-l 

Knrm,Nix,y) = ^ ipk{x)ipk{y), 

k=0 

which is then just the projection onto spanjc^Oj • • • ,^N-i}- Then the following formula 
holds: 

[-e^' ''^{I - Kurm,N)ix,y) ifs<t. 
Now introduce the rescaled process 

Af (t) = V2N^/''{xf{N-^/H) - Vm). 

Changing variables x i— t- ^^i/c, x + V'^N, y i— )• ^^i/a V + V^2iV in the kernel accordingly, 
a calculation gives 



Af (tl) < XI, . . . , Af (i„) < X„j = det(/ - fi^frmy)i2„,^,...,,„},ij) 

with _ 

.-e(^-*)^'^(/-AHrm,A)(2;,y) if s<t, 
where -^Hrm,Af(a;, y) = -^Hrm.A ( ^^1/6 + ^^^iV, ^/j^^ + \/2]V^ and the operator 

Ht^ = -A + X + 

The above rescaling corresponds to focusing in on the top curves of the Dyson Brownian 
motion. It is known that in the limit — )• 00, ^^Hrm.A converges to the Airy kernel Kp^i, 
while it is clear that converges to the Airy Hamiltonian H: 

(1.38) iJ = -A + x 

(i.e. H f{x) = —f"{x) + xf{x)). Putting aside precise convergence issues, the result is that 

~ , (e^'-^'^"KAi{x,y) ifs>t, 

(1.39) Ihn Klif^^^{s,x;t,y) = \ 

[-e'^'' ^'"{I - KAi)ix,y) ifs<t. 

The obvious question at this point is what is the relationship between this limit and the 
extended Airy kernel (1.21). It turns out that they are the same. This is a consequence of 
the following remark, which implies that the nice structure we saw in -fCnrm.Af survives in 
the limit: 

Remark 1.8. The shifted Airy functions (f)x{x) = Ai(2; — A) are the generalized eigenfunc- 
tions of the Airy Hamiltonian, as H(j)\ = X(px (we say generalized because (px ^ L'^(]R)). 
The Airy kernel i^Ai is the projection of H onto its negative generalized eigenspace. 
This is seen by observing that if we define the operator A to be the Airy transform, 
Af{x) := dz Ai{x - z)f{z), then Km = APqA*, where Pof{x) = l.<o/(x). 

In particular, e*^ is defined spectrally. Formally, its integral kernel is given by e*^ {x, y) = 
j^oo dXe~^'^ Ai(a; + A) Ai(y + A). The integral converges when t < by the decay properties 
of the Airy function, but it diverges when t > (it can be interpreted as 5x=y when t = 0). 
Nevertheless, in our formulas e*^ will always appear after Km when t > 0. This has the 
effect of restricting the integral to A > 0, which converges because the Airy function is 
bounded. 

As a consequence of the above discussion we obtain the following result, see [TW04]: 
Theorem 1.9. 

Af(i)T— ^^2(t) 

A— >oo 

in the sense of convergence of finite- dimensional ditributions. 
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The extended kernels which define the other Airy processes do not have exactly the same 
structure. One reason behind this is that, apart from A2 and ^1, the other processes are 
not stationary. The other reason is that, as we mentioned, in some cases (for example 
Airyi) the processes are obtained after further limiting procedures, which destroy part of 
the structure. Nevertheless, as we will see later, enough of this structure remains for our 
purposes. 

We have wandered a bit far from the main subject of this survey in the hope that the 
reader will get a feeling about why extended kernels appear naturally for our processes. In 
the rest of this article we will deal with formulas which are given as Fredholm determinants 
of certain operators acting on L^(M), as opposed L^({ti, . . . , t„} x M). In light of the above 
discussion, it is slightly surprising that such formulas should exist. 



2. Fredholm Determinants 
If K is an integral operator acting on H = {X, dfi) through its kernel 



(2.11 



iKf)ix) 



K{x,y)f{y)d^i{y), 



X 



we define the Fredholm determinant by 

00 



(2.21 



det(/ + Ai^) = 1 + 



A" 



n=l 



X 



X 



det [K{xi, Xj)]" . dfi{xi) ■ ■ ■ d/i(x„) 



If |i^(2;,y)| < B for all x,y, and /U is a finite measure, the Fredholm series (2.2) converges 
by Hadamard's inequality, 

|det(Ci,...,C„)| < ||Ci||---||C„|| 

where ||Ci|| denotes the Euclidean length of the column vector Cj, since the length of the 
column vector in [K{xi,Xj)]^j^^ is bounded by Bv}^"^, and hence the n-th summand in 
(2.2) is bounded by ^^"n'^/^. 

If one is not familiar with the definition (2.2) one might even wonder what it has to 



do with determinants. Take a matrix K 



d < 00, and consider the d x d 

determinant det(/ + XK). Clearly it is a polynomial of degree d in A, Yln=o ^nX"", and its 
coefficients are given by the rule 
for differentiating determinants. 



1 an 



det(/ + AK) |a=o- To compute this, use the rule 



dx det(Ci, ...,Cd) = Y^ det(Ci, . . . , dxCn, . . • , Q) 



n=l 



and the fact that, in our particular case, C„(A) = + XK.^n is linear in A and Cn(0) = en, 
the n-th unit vector. The result is 



det(I + XK) = 1 + X 



l<i<d 



l<i<j<d 



33 



+ A- 



l<i<j<k<d 



det 



K K 

^^3^ ^^33 
Kki Kk 



■3 



Kik 
Kjk 
Kkk 



+ ■■■ + X'^detK. 



Replacing the ordered sums with unordered sums gives a factor 1/n!, and setting A = 1 we 
can see that this is a special case of (2.2). Von Koch's idea [Koc92] was that this formula for 
the determinant was the natural one to extend to d = 00. Fredholm replaced the integral 
operator (2.1) on L^([0, l],dx) by its discretization [^K{^, to obtain (2.2), which 

he then used to characterize the solvability of the integral equation (/ + K)u = f via the 
non-vanishing of the determinant I + K. 
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One can of course imagine other, more intuitive definitions of the determinant. Perhaps 

(2.3) det(/ + K) = + A„) 

n 

where A„ are the eigenvalues of K, counted with multiplicity. Or 

(2.4) det(I + \K) = et=^i°g(i+^^) 
with the trace 

(2.5) tiK = J dn{x)K{x,x). 

Of course, these definitions require some smallness condition on K, but at least they 
make apparent the important fact that the determinant is invariant under conjugation 
det(/ + M-^KM) = det(I + K), or 

(2.6) det(/ + K1K2) = det(/ + K2K1), 

(usually referred to as the cyclic property of determinants) as well as the formula 

(2.7) dp det(/ + K(/3)) = det(I + K{I3)) tr((/ + K{P))-^dpK{P)) 
for K{I3) depending smoothly on a parameter /3. 

A more modern way to write (2.2) is 

00 

(2.8) det(/ + Ai^) = ^A"trA"(K) 

n=0 

where N^{K) denotes the action of the tensor product A® ■ ■ ■ ® A on the antisymmetric 
subspace oi H ® ■ ■ ■ ® H. If denotes the projection onto that subspace and C„ = 
P„A"(K)P„ then 

CnUl «> • • • ® /n) = ^ Sgn(c7)A/^(l) Af„^^) 

= ^ X] ^g'^(^) j' ■ 'j ^f^iVl) • • • dfJ^iVn) K{xi,y^(^i)) ■ ■ ■ K{Xn, ya(n))h{yi) ' ' ' fniVn) 

cr&Sn 

which shows that Cn is an integral operator with kernel det [K{xi, Xj)]"^.^^ and hence (2.8) 
is just a slick way to write (2.2). The advantage of (2.2) is that it can be used directly to 
define the Fredholm determinant for operators on a general separable Hilbert space, but 
we will not need this point of view here (see [Sim05] for more details). 

The natural notion of smallness for Fredholm determinants turns out to be the trace 
norm on operators 

||K||i :=tr\K\, 



where \K\ = \/K*K is the unique positive square root of the operator K*K. A (necessarily 
compact) operator with finite trace norm is called trace class. Using the Parseval relation, 
one can check that for such operators the trace can be defined as 

00 

n=l 

as it is basis independent. This works for operators on any separable Hilbert space, and in 
our setting it can be shown that this definition of trace coincides with (2.5) for K of trace 
class. The Hilbert- Schmidt norm \\K\\2 = y^tT{\K\^ is easier to compute, 

1/2 



|i^||2= dxdy\K{x,y)\^'^ 



and the relation between these norms and the more common operator norm HA'llop is 

lli^llop < lli^lb < lli^lll, 
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as well as 

||i^ii^2||i < ||i^i||2||i^2||2, pi^lli < Plloplli^lli, and Pi^lla < Plloplli^lb, 

all of which can be checked easily. Of course, in the latter two A need not be compact. 
The reason the trace norm is so useful is 

Lemma 2.1. 

1. (Lidskii's Theorem) If K is trace class then tr K = A„, where A„ are the eigen- 
values of K. It follows that the three definitions (2.2), (2.3) and (2.4) are equivalent. 

2. A det(/ + A) is continuous in trace norm. Explicitly, 

(2.9) |det(/ + E:i)-det(/ + i^2)| < ||i^i - /^alli exp(||iri||i + Hi^alli + 1). 

Lidskii's theorem is non-trivial and its proofs use heavy function theory, but (2.9) can 
be explained easily. Let f{z) = det(/ + ^{Ki + K2) + z{Ki — K2)), so that the left 
hand side of (2.9) is |/(|) — f{—^)\ < sup_x/2<t<i/2 l/'(OI- Cauchy's integral formula 
/'(^) = 2^^ f H^^' shows that sup_i/2<t<i/2 l/'(*)l < ji suP|^|<jj+i 1/(^)1- ^he eigenval- 
ues of A"(_fC) are Aj^ • • • Ai„, < • • • < in, so tr A"(i^) = X]ji<-<i„ ■ ■ ■ ^in and hence 
|trA'^(K)| < ^lli^lli, which implies 

|det(I + Air)| < e^ll^ll^ 
Therefore sup, .^^ \f{z)\ < exp{^\\Ki + K2II1 + {R + l)\\Ki - K2\\i) and taking R = 

I ^ I ^-H-H- 2 

||Ki -Ksllr^ gives (2.9). 



Examples. 

1. (Gaussian distribution) A trivial example is K{x,y) = ^7^=^- The operator is rank 



one, so if P, is the orthogonal projection from L (M) — )• L (s, 00) then by (2.4) we have 



det(/ - PsKPs)l^r,cI.) = 1 - tr P.i^P, = / 



dx 



ll-Kt 

Of course, the Gaussian here could be replaced by an arbitrary density. 

2. (GUE) Consider the Airy kernel Kxi{x, y) = dt Ai{x + t) Ai(y + t) and let Ait(x) = 
Ai(x + t) and H = —d"^ + x. Then -ffAi^ = — iAij, the Ait, t € M are generalized eigen- 
functions of and Kx\ is the orthogonal projection onto the negative eigenspace of 
H (see Remark L8). Using Ai"(x) = xAi(x), we have g /i(^+f)Ai'fa+t^-Ai'(x+f)Ai(^+t) ^ 

Ai(x + t)K\{y + t), which yields the Christoffel-Darboux formula 

, , A^x + t)A\{y + t)- A\{x + t)A]!{y + t) 

KAi[x,y) = . 

y - X 

To show PgKxiPs is trace class, write K^i = BqPqBq where 

(2.10) Bo{x,y)=Aiix + y). 
Then use 117^1^2111 < ||2||^2||2 to get 

/•OO POO 

(2.11) WPsKmPsWi < WPsBoPoWl < / Ai\x + y)dxdy, 

Jo J s 

which is finite by the following well-known estimates for the Airy function (see (10.4.59-60) 
in [AS64]): 

(2.12) |Ai(x)| < Ce^i^''" for x > 0, |Ai(x)| < C for x < 0. 
The GUE Tracy- Widom distribution is given by 

^gue(s) = det(/ - PsKAiPs)LHR,dx)- 
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On the face of it, it is not so obvious why such an expression would define a probabihty 
distribution function. From (2.11) it is clear that linis^oo det(I — PsKaiPs) = 1- Since 
PsKxiPs is a composition of projections, its eigenvalues satisfy 1 > Ai(s) > A2(s) > • • • > 0. 
Recall the min-max characterization of eigenvalues 

. , . . {f,PsKMPsf) 
Xk{s) = max mm — — , 

dimU=kfGU ifj) 

from which it is apparent that Xi{s) is non-decreasing as s decreases, and hence det(/ — 
PsK^iPs) = Y\i{l — Aj(s)) is non-increasing with decreasing s. In fact, Ai(s) 1 as 
s \ — oo since if / is in the negative eigenspace of H, {Psf, Kx\Psf) — (/, Kj^if) = (/, /). 
This shows that det(/ — PsKaiPs) \ as s \ — oo. 

3. (GOE) Fgoe(s) = det(/ - P,5o^s)L2(M,dx) where 5o(x,2/) is as in (2.10). The key to 
show that Bq is trace class in this case is the identity 

/oo 
dx Ai(a + x) Ai(6 - x) = 2~^/^ k\{2-^'^{a + h)) 
-oo 

(see, for example, (3.108) in [VSIO]). One defines Gi(x,z) = 2^/''^ Ai{2^/''^x + z)e'' and 
G2{z,y) = e"^ Ai(2V3y _ z) and notes that PsBqPs = {PsGi){G2Ps). Then (2.12) allows 
to show that each of the last two factors has finite Hilbert-Schmidt norm, yielding that 
PsBqPs is trace class. 

4. (Airyi process) Recall the Fredholm determinant formula (1.22) for the finite-dimensional 
distributions of the Airyi process. It turns out that the kernel iKf^^i inside the determi- 
nant is not trace class, basically because the heat kernel is not even Hilbert-Schmidt on 
^■^([5,00)) for s € M. Nevertheless, the series (2.2) defining the Fredholm determinant is 
finite in this case, because one can conjugate the kernel iKf^^l to something which can be 
proved to be trace class (see [BFP07]). 

The situation in the last example, where the natural expression for a kernel defines an 
operator which is not trace class, but which is conjugate to a trace class operator, arises 
often. Here by conjugacy we mean the following: two operators K and K are conjugate if 
there exists some invertible linear mapping U acting on measurable functions on X such 
that K = UKU~^. Observe that such a pair of operators have the same Fredholm series 
expansion (2.2), i.e. det{I + K) = det{I + K). This allows to extend the manipulations on 
Fredholm determinants to operators which are conjugate to trace class operators, provided 
that one is careful in keeping track of the needed conjugations. 

The reason we start with the Fredholm expansion (2.2) is that this is the way the 
determinant usually arises from combinatorial expressions. Sometimes the kernels are not 
trace class, but this should not bother us so much as long as some version of the formal 
expression can be shown to converge, for instance as in Example 4 above. Often, it is 
genuinely difficult to show that the resulting expressions define a probability distribution, 
and we only know it because they arose this way. 

3. Boundary Value Kernels and Continuum Statistics of Airy Processes 

3.1. Boundary value kernel formulas for finite-dimensional distributions. Recall 
the formula (1-19) for the finite-dimensional distributions of the Airy2 process. It is given 
in terms of the Fredholm determinant of what we call an extended kernel, that is, (the 
kernel of) an operator acting on the "extended space" L'^iih, ...,tn}xR). Although 
such formulas have been very useful in the study of models in the KPZ class, they suffer 
from two problems. First, if one wants to take the number n of times ti to infinity, a 
big difficulty appears in the fact that these formulas involve Fredholm determinants on 
the Hilbert space L^{{ti, . . . ,tn} x M), and thus the space itself is changing as n grows. 
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Second, these formulas are useful for computing long range properties of the processes (for 
instance an asymptotic expansion of the covariance of A2{s) and A2{t) as \t — s\ — )• oo, see 
[Wid04] ) , but are not suitable for studying short range properties such as regularity of the 
sample paths. 

The second type of Fredholm determinant formula, which is the one we will use for 
most of the rest of this article, was actually introduced as the original definition of the 
Airy2 process by Prahofer and Spohn [PS02]. It is given as follows: for ti < • • • < t„ and 

Xi, . . . ,Xn €R, 

(3.1) F{A2{tl) <Xi,...,A2{tn) <Xn) 

= det(/-KAi + P..e(*-*^)^P..,e(*-*^)^...P.„e(*"-*i)^KA^^ , 

where K^i is Airy kernel (1.9), H is the Airy Hamiltonian (1.38) and Pa denotes the 
projection onto the interval (— oo,a]: 

Note that the Fredholm determinant is now computed on the Hilbert space 

L2(m) instead 

of L^({ti, . . . ,t„} X M), which makes taking n — )• cxd at least feasible. Note also that the 
time increments ti — tj+i appear explicitly in the formula, which explains why this formula 
will be more suitable for the study of short range properties. Another advantage of this 
formula is that it makes apparent that A2 is a stationary process. 

The equivalence of (1-19) and (3.1) was derived formally in [PS02] and [PSll]. The proof 
in [PSll] is based in the following idea. As we explained in Section 1.5.1, the extended 
kernel formula (1.19) can be thought of as the determinant of an n x n matrix whose entries 
are operators acting on (M) . By rewriting this operator as a sum of an upper-triangular 
part and lower-triangular part and using algebraic properties of the determinant and the 
algebraic relationships between the different entries of this matrix, [PS02] showed that the 
determinant equals the determinant of an operator- valued matrix I + G such that only the 
first column of G is non-zero. Therefore det(/ -|- G)i2({j^_ j^}xr) = det(/ -|- G'i^i)^2(]g) (to 
see this simply pretend the operators in the determinants are matrices), and an explicit 
calculation of Gi^i yields (3.1). 

The argument given in [PSll] which we just sketched is almost a complete proof. There 
are nevertheless some subtleties. For example, it is not a priori obvious that for s,t > 0, 
e~^^ can be applied to the image of Pae~^^ . Moreover, in order to manipulate Fredholm 
determinants one needs to check that certain analytical conditions are satisfied (see Section 
2). The technical details are discussed in [QR12a], which in fact shows that a formula 
analogous to (3.1) holds for the Airyi process as well. It is given as follows: for < • • • < i„ 
and xi, . . . , x„, G M, 

(3.2) P(^i(ti) <xi,...,^i(t„) <x„) 

where Bq is given by the kernel Bq{x, y) = Ai(x-l-y) defined in (1.13) and A is the Laplacian 
operator. Observe that in all but the last factor of the form e*^ in the above formula it 
holds that s > 0, in which case e''^ is the usual heat kernel. This kernel is ill-defined for 
s < 0, but it turns out that in this case the operator c'^'^Bq makes sense if defined via the 
integral kernel 

(3.3) e'^Boix, y) = e2«V3+^(x+y) + y + ^2)^ 

What we mean by this is that if s, t > then, with this definition the semigroup property 
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KT\s,x;t,y) 



As we will see in Section 3.2, it is fruitful to think of the operator appearing in (3.1) 
as the solution of certain boundary value problem, so we will refer to formulas like this as 
boundary value kernel formulas. By using (3.3) one can rewrite the definition (1.23) of the 
extended kernel for Ai as 

'e(*-'^)^i?o(x,y) ifs>t, 

-e(*-^)^(/-i?o)(x,y) \is<t. 

It becomes clear then that both the extended kernel formula and the boundary value kernel 
formula for Airyi are obtained from the corresponding formulas for Airy2 by substituting 
H with —A and Kp^\ with Bq. It turns out, as shown in [BCR12a], that the necessary 
structure behind these formulas hold for a much wider class of processes, including for 
instance the stationary GUE Dyson Brownian motion and non-stationary processes like 
the Airy2_!.i process, and the Pearcey process [TW06]. For example, for Airy2_s.i one has 
[BCR12a] 

(3.4) P(^2^i(ti) < XI, . . . ,^2^l(t„,) < Xn) 

— f\p.i(T—K*^ -4- P- p(*2-<i)Ap_ (t„-t„_i)Ap_ {ti-t„)Azfti \ 

where Xj = Xj — tflt^Ko and K2^]^(x,y) = i^2^]^(t, x; t, y) with as in (1.24). 

Interestingly, it is shown in [BCR12a] that in a setting corresponding to discrete non- 
intersecting paths, analogous boundary value kernel formulas can be obtained directly 
from applying the Karlin-McGregor formula (1-37) (or rather its combinatorial analogue, 
the Lindstrom-Gessel-Viennot Theorem [Lin73; GV85]), bypassing the direct application of 
the Eynard-Mehta Theorem. In the case of the Airy2 process, a suitable limit of a discrete 
family of non-intersecting should lead to (3.1) (cf. (1.39)). Such a procedure does not 
seem to work for the Airyi process. In fact, in that case the determinantal process used 
to derive (1.23) is signed (in the sense that the measure defined by the analog of (1.35) is 
signed), see [BFPS07], and hence it is not clear how to associate directly to it a family of 
non-intersecting paths. 

As we will see below, the integral kernels of the operators appearing inside the Fredholm 
determinants in (3.1), (3.2) and (3.4) can be expressed simply in terms of hitting probabili- 
ties of Browian motion. In other words, hitting probabilities of curves by A2, Ai and A2^i 
can be expressed in terms of Fredholm determinants of the analogous hitting probabilities 
for Brownian motion. Given the above discussion (and the discussion in Section 1.5.2), this 
is not entirely surprising in the case of A2, as it follows from the non-intersecting nature 
of systems of Brownian paths that can be used to approximate A2- For the same reason, 
it is surprising in a sense that the same structure is present in ^1. 

3.2. Continuum statistics and boundary value problems. Consider the following 
problem: compute the probability that, inside a finite interval r], the Airy2 process lies 
below a given function g. The obvious way to proceed is to take a fine mesh i = ti < t2 < 
■ ■ ■ < tn = r of the interval [i, r], take Xj = g{ti), and attempt to take a limit as n — ?• c« in 
the formula for the finite dimensional distributions of A2, 

(3.5) F{A2{h) < g{ti), . . .,A2{tn) < g{tn)) 

= det (l - i^Ai + ^'.(0^^*^'*^^''^.(*2)e(*^-*^)^ • • • P,(t„)e(*"-*^)^i^Ai) . 

Here the Fredholm determinant is computed on the Hilbert space -L^(M), which we will 
omit from the subscript in the sequel. By Theorem 1.1, A2 has a continuous version, 
and hence lim^^oo P(^2(ti) < 5(^1), • • ■ ,^2(tn) < 9{tn)) = HMt) < 9{t) for t G %r]). 
To study the right hand side of (3.5) we need to compute the limit of the operator appear- 
ing inside the determinant. Observe that the last exponential equals e^^~^^^ , and hence 
does not depend on n. On the other hand, for s < t the operator e^*"*)-'^ can be thought 
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of as mapping a function / to the solution u(t, •) at time t of the PDE dfU + Hu = with 
initial condition u{s, ■) = /(•). Therefore the operator 

can be thought of as solving the same PDE (backwards in time) on the interval [£, r] with 
the additional condition that all the mass above g{ti) is removed at each of the discrete 
times tj. Note that the PDE is solved backwards because, if we apply this operator to a 
function on its right, we first apply Pg(t„)j then g(tn-i-t„)H ^ then Pg^tn-i)-! ^^^^ so on. Since 
solving the PDE dtu + Hu = forward or backwards in time gives the same answer, if 
we want to think of (3.6) as being solved forward in time, all we need to do is reverse the 
order in which the g{ti) appear. The result is the following. Given g G ff^([£, r]) (i.e. both 
g and its derivative are in Lp'{[l,r\)), define an operator 0^£^] acting on L^(]R) as follows: 
0|^^j/(-) = u{r, •), where u(r, •) is the solution at time r of the boundary value problem 

dtu + Hu = for X < g{t), t e {i,r) 

u{e,x) = f{x)i^^g(^(,) 

u{t, x) = for X > g{t). 
Further, define g{t) = g(i + r — t). Then 

Since the convergence holds in trace class norm, (3.7) can be used to answer the question 
with which we started this subsection: 

Theorem 3.1 ([CQR12], Theorem 2). 

(3.8) F{A2{t) < g{t) for t G [£,r]) = det(l - Km + e^./'-'^'^KM 



Observe that we have written g instead of g in (3.8). We may do this because the Airy2 
is invariant under time reversal, so we can replace g hy g on the left hand side. 

The limit (3.7) is proved in Proposition 3.2 of [CQR12] (in fact only along the dyadic 
sequence = 2^, but this is enough for deducing Theorem 3.1). The proof is based on the 
following probabilistic representation of the solutions of the above boundary value problem: 
if 0^^j,j(x,y) denotes the integral kernel of 0^£^], then 



(3.9) ef,,.j(x,y) 



9 „^ ^ ^ex-ry+ir^-e'^)/3 



g-(a;-y)2/4(r-f 



where the probability is computed with respect to a Brownian bridge h{s) from x — l"^ at 
time £ to y — at time r and with diffusion coefficient 2. This formula is Theorem 3 of 
[CQR12] , its proof is based on an application of the Feynman-Kac and Cameron-Martin- 
Girsanov formulas. 

The argument that proves Theorem 3.1 can be adapted to obtain a similar result for 
Airyi. Fix i < r. Given g G H^{[i,r]), define an operator A^^^j acting on L^(M) as follows: 

A^^j/(-) = u{r, •), where u{r, ■) is the solution at time r of the boundary value problem 

dtu - Au = for X < g{t), t G (^, r) 
u{£,x) = f{x)l^^g(_e) 
u{t,x) = for X > g{t). 
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Theorem 3.2 ([QR12a], Theorem 4). 

(3.10) ¥{Ai{t) < g{t) for t G [e,r]) = det(/ - Bq + Af.^^je-^^-^^^^o) • 

Although the proof of this result is similar to the proof for the Airy2 case, the argument 
is a bit more involved because, as written, the operator in the determinant is not trace 
class, so one needs to conjugate appropriately. Of course, similar arguments should allow 
one to obtain continuum statistics formulas for other processes for which boundary value 
kernel formulas are available (see [BCR12a] for the case of stationary GUE Dyson Brownian 
motion) . 

The operator A^^ also has a simple representation in terms of Brownian motion (see 
[QR12a]), which has recently been used in [FF12] to verify numerically the experimental 
values obtained in [TS12] for the persistence probabilities of Airyi. The negative persistence 
exponent is defined by 

¥{Ai{t) <m, < t < L) ~ e""--^ 

where m is the mean of -Fgoe- Takeuchi has measured k_ ~ 3.2 it 0.2 in computer 
simulations of the Eden model [TS12]. Ferrari and Frings [FF12] have computed numerically 
(3.10) finding 

K w 2.9, 

which is fairly close. Note that Takeuchi has also measured the positive persistence prob- 
abilities F{Ai{t) > m, < i < L) ~ g-i^+L interesting question is whether there exists 
a simple enough mathematical formula to check such a thing. 

4. Applications 

In this section we will describe some applications of the boundary value kernel formulas 
for Airy processes which were introduced in the previous section. The first two applications 
refer to asymptotic statistics for directed polymers and LPP, while the next two involve 
respectively the Airyi and Airy2-!.i processes. 

4.1. Point-to-line LPP and GOE. Recall the variational formula (1.25) relating the 
Airy2 process with the Tracy- Widom GOE distribution: 

(4.1) P( sup {^2(2;) - x^} < m ) = FGOE(4^/^m). 

As we explained in Section 1.2.3, Johansson's proof [Joh03] was very indirect, relying on 
the convergence of the spatial fluctuations of point-to-point LPP to A2 together with (1.15) 
and (1.17). 

A direct proof of this variational formula was provided in [CQR12], based on the contin- 
uum statistics formula given in Theorem 3.1. An interesting consequence of this derivation 
was that it allowed to identify the factor of 4}^^ on the right hand side of the identity, 
which had been lost in Johansson's argument in the process of translating between the 
available results at the time (see Section 2 of [CQR12] for an account of how to get the 
correct factor directly from LPP). 

We will explain next the derivation of the formula, skipping some details. We rewrite 
the desired probability as 

lim ¥(A2{t) <m + t^\/te [-L,L]) . 

L— s>oo 

For fixed L > 0, Theorem 3.1 implies that this probability is given by 

(4.2) det(/ - Km + e^e^^^i^Ai) , 
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where 

The nice thing is that the choice of g{t) = + m is the simplest possible from the point 
of view of explicit calculations, because it cancels exactly the parabola appearing on the 
right hand side of (3.9). The probability appearing in that formula is then reduced to the 
probability of a Brownian bridge staying below level m, and this is easy to compute using 
the reflection principle (method of images): 



^h{~L)=x-L^h{L)=y-L^\;'^^) < ^ OH [-L,L] 

(4.3) = 1 _ Pg(_^)^^_i2,S(i)=j^_i2 (b{s) > m for some s £ [-L, L] 

= 1 - € 



'{x-m-L'^)(y-m-L'^)/2L 



(we leave the simple computation to the reader, alternatively see page 67 of [BS02]). 
Putting this back in 0^ gives 

where Rl is the reflection term 

(4.5) RUx, y) = _^e-("+J^-2™-2i')V8i-(-+?/)^+2iV3. 

The e^^^^ in the first term in comes from the 1 in (4.3) and appears by either reversing 
the use of the Cameron-Martin-Girsanov and Feynman-Kac formulas in the derivation of 
(4.4) or by an explicit computation of the integral kernel of e~^^~^^^ as 

— (x—'ti)'^ /4(r—£) 



V^4^(r - e) 

Referring to (4.2), we have by the cyclic property of determinants (2.6) and the identity 
e^^^Kxi = (e^^^Ai)^ (which follows from Remark 1.8) that 

(4.6) P(^2(i) <t^ + mioTt£ L]) = det(/ - + e^" K^iQie^" K^i) ■ 

To obtain the L — t- oo asymptotics, we decompose Ql so as to expose the two limiting 
terms, as well as a remainder term VI l: 

(4.7) QL = e-^^^ -Rl + ^l, 

where Q.l = {Rl - Pm+L'^RLPm+L^) - (e"^^^ - Pm+i'^e''^^^ Pm+L-^) ■ It is shown in 
[CQR12] that 

(4.8) lle^^iTAiJ^Le^^KAilL > 0. 

The proof amounts essentially to asymptotic analysis involving the Airy function. In view 
of this fact and the decomposition (4.7), and since e^^ K Xie~^^^ e^^ K x\ = Kxi, we see 
that the key point is the limiting behaviour in L of 

e^^i^AiiiLe^^^Ai. 

To explain how this last product can be computed we will proceed in a slightly formal 
manner through an argument based on the Baker-Campbell-Hausdorff formula, as done for 
a related problem in [QR12b] (see Section 4.4). Since i^Ai is a projection and H leaves 
invariant, we will pretend that e^^ and Kx\ commute, so we have to compute the limit of 
e^^ Rlg^^ . Define the reflection operator Qm by 

Qrnfix) = /(2m - x). 

Then the operator Rl defined in (4.5) can be rewritten as 

(4.9) Rl = e(2^')/3e-^«£>^+i2e2^^e-^« = e(2^')/3e-^«e^^^.^+i2e^^e-^«. 
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Here e^^ stands for a generic variable) denotes the multiplication operator (e^^/)(x) = 
e^^f{x). The second equality follows from the reflection principle applied to the heat kernel. 

The following identities will be useful, where [•, •] denotes commutator: 

[H, A] = A] = -2V, [H, V] = [C, V] = -/, [H, ^] = -2V. 

If A and B are two operators such that [A, [A,B]] = cil and [B, [A,B]\ = C2I for some 
ci, C2 G M, then the Baker-Campbell-Hausdorff formula reads 

gAgS ^ ^a+B+\[A,B] + ^[A,[A,B]]-UB,{A,B]]_ 

Using this we have 

Using the Baker-Campbell-Hausdorff formula again we deduce that 
while an analogous computation yields 

gLAg-L€gL// ^g-L3/3g-L2V^ 

Employing these identities on the right hand side of (4.9) yields 

Since e*"^ is the shift operator (e^^/)(a;) = f{x + r), we have e^^ Qm = gme~^'^ = Qm-r/2^ 
and we obtain 

Remarkably, the result does not depend on L. The conclusion from using this, (4.8) and 
(4.6) in (4.7) and taking L — )• 00 is that 

(4.10) P(A(t) < + m for all t G M) = det(/ - KAiQmKAi) ■ 

The use of the Baker-Campbell-Formula in the derivation of this identity can be replaced 
by an explicit integral calculation (see the proof of Proposition 1.3 of [CQR12]). 

To finish our proof of (4.1) we need to show that the right hand side of (4.10) equals 
i^GOE(4^'^^m). Recall the definition of the kernel B(){x,y) = A\{x + y) and observe that 
i^Ai = BqPqBq. Recall also that the shifted Airy functions form a generalized orthonormal 
basis of L^(M) (see Remark 1.8), which implies that Bq = I. Therefore we can use the 
cyclic property of determinants (2.6) to deduce that 

det(/ - KAiQmKAi) = det(/ - PoBoQmBoPo) . 

Now 

/oo 
dX Ai{x + A) Ai(2m - \ + y), 

and using the identity (2.13) we deduce that 

(4.11) BoQmBoix, y) = Bm{x, y) := 2'^/^ Aii2-^/\x + y + 2m), 
and thus 

det(/ - KAiQmKAi) = det(/ - PoBmPo) ■ 

Performing the change of variables x 1— )• 2^/^x, y 1— >• 2^/^y in the series defining the last 
Fredholm determinant shows that the determinant on the right hand side of (4.10) equals 
det(I- Po54i/3^Po), which is FGOE(4i/='m) by (1.14). 
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4.2. Endpoint distribution of directed polymers. In the setting of geometric LPP 
(see Section 1.2.3), consider the random variable 

KN = mm{ k G {-N, ...,N}: sup = sup > . 

^ j=-N,...,k j=-N,...,N j 

kn corresponds to the location of the endpoint of the maximizing path in point-to-line 
LPP. 

Interest in the scaling properties and distribution of this random variable goes back at 
least to the early 1990's. One can also consider the analogous random variable in the setting 
of directed random polymers, but due to the KPZ universality conjecture one expects that 
the asymptotic behavior and statistics are the same as in LPP. Mezard and Parisi [MP92] 
considered the polymer case and derived non-rigorously the scaling relation 

\kn\ ~ iV^/^ 

(c.f. (1.6)). In view of this we define the rescaled endpoint 

where cs is the constant appearing in (1.18). Recalling the definition of the rescaled point- 
to-point last passage time (1.18) as the linear interpolation of the values given by 

FP™''*(t) = — 3_ rLP™°t(iV + c^^N~^/% N - c^^N-'^/h) - ciN 

for t such that c^^N^'^/^t G {-iV, . . . , iV} we deduce that 

Tn = min< t G M: sup -fr^°™*(s) = supif^°™*(s) 
I s<t seR 

Recalling that H]\f(t) converges to A2{t) - by Theorem 1.1 it becomes clear that T/v 
should converge to the point where A2{t) — attains it maximum. In fact, this is what 
Johansson proved, although he had to make a (very reasonable) technical assumption on 
the Airy2 process which he was not able to prove with the tools available at the time: 

Theorem 4.1 ([Joh03]). Assume that the process A2{t) — t^ attains its maximun at a 
unique point and let 

T = argmax {A2{t) — t^}. 

tm 

Then 

^^Tr 

iV-5>oo 

in the sense of convergence in distribution. 

Although the result is of course very interesting, as it shows that the limiting endpoint 
distribution exists (under the technical assumption) , it gives no information on the distribu- 
tion of 7". Quoting Johansson [Joh06], for all we know T could be Gaussian. Nevertheless, 
from KPZ universality one expects that this is not the case. For example, Halpin-Healy 
and Zhang [HHZ95] conjectured on the basis of analogy with the argmax of a Brownian 
motion minus a parabola (for which one has a complete analytical solution, see [Gro89]), 
that the tails of T decay like e"*^* , which of course rules out Gaussian behavior. 

It turns out that the distribution of T can be computed explicitly through an argument 
based on the continuum statistics formula of Theorem 3.1. This was done in [MQR12], 
where in fact the joint density of 

T= argmax {^2(t) - i^} and M = m.a-^ {A2{t) - t^] 
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was computed. Moreover, the argument implies that the maximum of A2{t) —t^ is attained 
at a unique point, thus completing the proof of Theorem 4.1. The uniqueness of the max- 
imum was also proved slightly earlier by Corwin and Hammond [CH12] using completely 
different techniques, and a proof for general stationary processes is now available [Piml2]. 

The computation is as follows. For simplicity we will assume the uniqueness of the 
maximizing point of A2{t) — t^, and will explain later how the uniqueness can actually 
be obtained from this argument. Let {Ml^Tl) denote the maximum and the location of 
the maximum of A2{t) — restricted to i S [— L,L], and let be the joint density of 
{M.l-,Tl)- By results of [CHll], the joint density f{m,t) of 7W, T is well approximated by 
fL{m,t), 



By definition, 



f{t,m) = lim fL{t,m). 
L—^oo 



fL{t,m) = lim lim ^¥{Ml G [m,m + e], Tl G [t,t + 6]), 
5^0 £-s>o eo 



provided that the limit exists. The main contribution in the above expression comes from 
paths entering the space-time box [t, t + 5] x [m, m + e] and staying below the level m 
outside the time interval [t, t + 5]. More precisely, if we denote by 5 and D^^s the sets 

De,s = {A2{s) - s"^ < m, s £ [t,t + 6f, A2{s) - s"^ < m + e, s£[t,t + 5], 

A2{s) — S [m, m + e] for some s S [t, t -|- (^]|, 

and 

Ds,s = 1^2(5) — s"^ < m + £, s € [— L, L], A2{s) — G [m, m + e] for some s G [t, t -|- 
then 

D,5 C {Ml €[m,m + e],TL £ [t,t + 6]} ^D,^s- 

Letting f{t,m) = lim5_^o li™£-s>o ^IF'(iZ£,5) and defining f{t,m) analogously (with 0,, ^ 
instead of s) '^^ deduce that f{t,m) < f{t,m) < f{t,m). In what follows we will 
compute f{t, m). It will be clear from the argument that for /(t, m) we get the same limit. 
The conclusion is that 

/,(t,„^) = hmlimlF(^,). 
We rewrite this last equation as 

fL(.t,m) = lim lim ^[P(^2(s) < K^s{s), s G [-L, L]) -F{A2{s) < ho^s), s G [-L,L])], 
where 

he,s{s) = s"^ + m + els(z[t^tj^s]- 

These two probabilities have explicit Fredholm determinant formulas by Theorem 3.1. We 
get, using the cyclic property of determinants as in (4.6), 



fr (t, m) = lim lim ^ 



det(l - Km + e^^i^AiOj':-^ .^je^^i^Ai 



det ( / - Km + e'-'^KMeili ne^^i^Ai 



The limit in e becomes a derivative 

(I-KM + e^^K^^e':^\.e^"K^^ 

13=0 



fUt, m) = lim ^ dp det (l - Km + e^^i^AiOjl^^^ ^^jC^^^Ai 
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• hm - tr 

(5-5-0 



(/ - Km + e^^i^Aie[_!;^^^]e^^KAi)-'e^^i^A 



(4.12) fL{t,m) = detil - Km + e'^'' KmQ7Jj ne^^i^Ai 



LH - 



^1„LH - 



/3=0 



e^^i^Ai 



One has to check here that the required limits hold in trace class norm, see [MQR12]. Note 
that /io,5 = Qm, where is the parabolic barrier 



s + m, 



so in particular the determinant and the first factor inside the trace do not depend on 6. 
We know moreover from the arguments in Section 4.1 that 

^lirn^ (/ - Km + e'^''KMQ\yL,Lf''"KA) =1- KMQmKM 
in trace norm. In particular, we have 

^lim^ det(l - Km + e^^KAiOj:!^ .^je^^KAi) = FGOEi^^'^'rn). 

The next step is to compute dj^Q^^'^ ^ 1/3=0- Recalling that h^ sis) = gm{s) = + m 
and also h^^sis) = gm+eis) for s £ [t,t + S] we have, by the semigroup property, 



e 



e 



ho.s 



-L,L] 



L,t\ 



C!t9m 



Computing the desired derivative involves just the middle bracket, which we note corre- 
sponds to the same boundary value problem as in Section 4.1, only at two different levels m 
and m + e. Since we have explicit formulas, the derivative can be computed explicitly. The 
computation is slightly tedious, and the only delicate part is to justify that the necessary 
limits occur in trace class norm, we refer to [MQR12] for the details. 

Going back to (4.12), we recall that the trace is linear and continuous under the trace 
class norm topology, so in view of the preceding discussion we have 



(4.13) lim fL{t,m) 
L-^oo 

= FGOE(4^/^m)tr 



(/ - KMQmKMr^ hm lim ^e^^i^Ai 



/3=0 



e^^i^Ai 



Once again we need to compute limits, again taking care that they hold in trace class norm 
as necessary. We skip the details and just write down the result, (2.12) in [MQR12]: 



1 



hm hm ^c^^Kai 

L-^oo <5— s>0 



L,L] 



/3=0 



e^'^KM 



where 
and 



V^t,„(x) = 2e*'+('"+^)* 



Ai\m + t'^ + x) + t Ai(m + + x) 



(we remark that we have written these formulas in a slightly different way compared to 
[MQR12], but the reader should have no problem translating between the formulas). The 
limit in 6 is relatively straightforward, while the limit in L involves an argument similar to 
the one used in Section 4.1. Using this formula in (4.13) the trace becomes 

tr [(/ - KAigmKAi)'^"^] = {{I - KAigmKAi)~^BoPo^t^rn, BoPoTp-t,m) , 

where (•, •)-^ denotes inner product in L^(M). 
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It only remains to simplify the expression. We first use (4.11) and the facts that Kji^\ 
'1 = 1 and Bl 



BqPqBq, Bq = I and Bq = Bq to write 



((/ - KAigmKAir^BoPoi't,m, Bo^'oV'-t.m) = ((^ " B^P^BraPM'^ B^P^i't,m, BoPoij^t,m) 

= {Bo{I - PoB^Po)-^PoiJt,m, BoPoi^-t,m) 
= ((/ - PoBmPo)-^Poi^t,m, PQ^-t,m) ■ 

Next we introduce the scaling operator Sf{x) = One can check easily that 

5-1 = and that Pq commutes with S and S~^. We also have SBmS~^ = -641/3^. 

Thus writing fh = 2~^/^m we get 

= {S-\I - PQB2rhPor^PoS^t,m^Poll>-t,m) 

which is equal to 2"^/^7(t, 

Using this formula in (4.13) yields the joint density of T and M. Define the resolvent 
kernel 

?m(a:,2/) = {I - PQBmPo)~^{x,y) 

and, for t, m G M, define 

■^t,m{x,V) = 2^/3V'i,m(2'/^2;)V-i,™(2l/3y) 

and 



00 fOO 



7(t,m) = 2V3/ dyV'-t4-i/3m(2'/'a;)w(x,y)V'i4_i/3„(2i/3y). 
JO Jo 

Theorem 4.2 ([MQR12], Theorem 2). The joint density f(t,m) ofT and M is given by 

/(t,m) =7(t,4i/='m)FGOE(4^/=^m) 
(4.14) , 

= det(/ - PoB^^/s^Po + Po^t,mPo) - PGOE(4^/3m). 

To see where the second equality in (4.14) comes from, observe that 'y(t,4}^^m) equals 
the trace of the operator {I — PoB^i/s^Pq)'^ PQ"^t,mPo and that "^t,m is a rank one operator. 
The identity now follows that from the general fact that for two operators A and B such 
that B is rank one, one has det(I - A + B) = det(/ - A)[l + tr ((/ - A)-^B)] . 

Integrating over m one obtains a formula for the probability density /end(^) of T. Un- 
fortunately, it does not appear that the resulting integral can be calculated explicitly, so 
the best formula one has is 



/end(i) = / dmf{t,m). 
J —00 

One can readily check nevertheless that /cnd(i) is symmetric in t. The second formula for 
f{t, m) is suitable for numerical computations, using the numerical scheme and Matlab 
toolbox developed by Bornemann in [BorlOa; BorlOb] for the computation of Fredholm 
determinants. Figure 3 shows a contour plot of the joint density of Ai and T, while Figure 
4 shows a plot of the marginal T density. 

As we mentioned, interest in this problem dates back at least two decades. In particular, 
there has been a resurgence of interest in the last couple of years. An alternative way 
to obtain the Airy2 process is as a limit in large N of the top path in a system of N 
non-intersecting random walks, or Brownian motions, the so called vicious walkers [Fis84] 
(this is of course related to the setting presented in Section 1.5.2). [SMCRF08; Fei09; 
RSIO; RSll] obtained various expressions for the joint distributions of M and T in such 
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Figure 3. Contour plot of the joint density of Ai and T. 



a system at finite A^. [FMSll] obtained the i*GOE distribution from large N asymptotics 
non-rigorously, and furthermore made connections between these problems and Yang-Mills 
theory. But for several years people were not able to perform asymptotic analysis on the 
formulas obtained for T at finite A^. 

After [MQR12] appeared, Schehr [Schl2] succeeded in extracting asymptotics from the 
vicious walkers formula, and obtained an alternative formula for f{t,m). His formula is 
given as follows. The Painleve II equation (l.lO)-(l.ll) has a Lax pair formulation 



(4.15) 



ds 



for a two-dimensional vector $ = <I>(C,s), where the 2x2 matrices A = A{(^,s) and 
B = B{C,, s) are given by 



-4(2 - s - 2g2 + 2q' 



4(2 + s + 2g2 + 2q' 
-Kq 



and 5(C) s) 



q C 
-C -q 



The compatibility of this overdetermined system implies that q{s) solves Painleve II. Now 



let $ = ( 1^ ) be the unique solution of (4.15) satisfying 



$i(C;s) = cos( ^c' + <) + o(r^^ 



as C — )• iboo for s € M. The formula of [Schl2] is 

(4.16) ^{t,m) = ^{\2/st,h_^2/3t)L2(m,oo) 

where 

l-OO 

ht{x)= / dCCMC,x)e-'^\ 



sin( ^c' + <) + o(r'^ 



Although Schehr's argument is non-rigourous, a later paper of Baik, Liechty, and Schehr 
[BLS12] proved directly the equivalence of the formula of [Schl2] and (4.14), thus estab- 
lishing the validity of (4.16) based on Theorem 4.2. 
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Figure 4. Plot of the density of T compared with a Gaussian density with 
the same variance 0.2409 (dashed Une). The excess kurtosis E(7'^) /E(T^)^ — 
3 is -0.2374. 



Before turning to the tail behavior of T, let us briefly explain how the uniqueness of 
the maximizer of A2{t) — can be established directly from the argument we described 
above. In the derivation of the formula we assumed that the maximum of A2{t) — 
is obtained at a unique point. However, it is not necessary to do this. In fact, if one 
follows the argument without this assumption, one ends up with a formula for what is in 
principle a super-probability density, i.e. a non-negative function /(t, m) on M x M with 
J^^^dmdt f{t,m) > 1, and in fact one can see from the argument that 

/ ,m it m m) = expected number of maxuna of Mt) - 

JRxR 

Recall that from (4.1) that the distribution of A4 is given by a scaled version of -Fqoe- A. 
non-trivial computation (see Section 3 of [MQR12]) gives 

/oo 
dt/(t,m) = 4i/3F^OE(4'/'"^)- 
-oo 

This shows that f{t,m) has total integral 1, which can only be true if the maximum is 
unique almost surely, since the global maximum is attained at at least one point. 

We mentioned earlier the conjecture that T should have tails which decay like e~^*^ (see 
e.g. [HHZ95]). This can be proved using the techniques described in this review: 

Theorem 4.3 ([QR12c; Schl2; BLS12; CHll]). There is a c> such that for every k> ^ 
and large enough t, 

[CHll] had obtained the e"*^*^ decay for some c > 0. The statement we included here is 
the one appearing in [QR12c]. In fact, Schehr's formula and its validation in [BLS12] later 

4 +3 

yielded a lower bound that matches the e~3 behavior of the upper bound, so we know 
now that | is the correct exponent. A precise asymptotic expansion of P(|7'| > t) based 
on that formula has recently been obtained in [BL12]. The reason why [QR12c] obtained a 
slightly worse lower bound is technical, and arises from the fact that the explicit formula 
(4.14) for f(t,m) is not useful for providing a lower bound, and instead one needs to use 
a different argument. On the other hand, the upper bound can be obtained directly from 
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(4.14). In fact, the second formula expresses this joint density as the difference of two 
Fredholm determinants, so we may use (2.9) to estimate the difference, and then all that 
remains is to show that this estimate can be integrated in m. See [QR12c] for more details. 

4.3. Local behavior of Airyi. As we mentioned, the boundary value kernel formulas 
introduced in Section 3.1 are better adapted than the standard extended kernel formulas 
to study short range properties of the processes. An interesting application is the following: 

Theorem 4.4 ([QR12a], Theorem 2). The Airy\ process Ai and the Airy2 process A2 
have versions with Holder continuous paths with exponent ^ — S for any 6 > 0. 

Continuity was known for A2 (see Theorem 1.1) but not for Ai. The Holder ^— con- 
tinuity for A2 also follows from the work of Corwin and Hammond [CHll]. Their proof 
is based on a certain Brownian Gibbs property for the Airy2 line ensemble (an infinite 
collection of continuous, non-intersecting paths, the top line of which is A2), and as such it 
cannot be extended to Airyi, given that no analog of the Airy2 line ensemble is known in 
the flat case. This regularity is expected to hold in fact for all the Airy processes in view 
of the fact that they are believed to look locally like a Brownian motion (see Section 1.3). 
Analogous results have recently become available for the solutions of the KPZ equation at 
finite times with certain initial conditions [Haill; QRll; CH12]. 

The proof of Theorem 4.4 is based on an application of a suitable version of the Kol- 
mogorov criterion. In the Airyi case, it involves studying a truncated version of the process, 
Af{t) = ^i(t)l|^^(i)|<M + Ml_4^(t)>M - Ml_4j(t)<_M and then proving the following es- 
timate: for fixed 5 > 0, there is a io G (0, 1) and an no G N such that for < t < to, n > no 
and M = (31og(t~(i+")))^''^ we have 

E([ylf(t)-^f(0)]'") <ct^+(i-'^)" 

where the constant c > is independent of 6, no and to. The proof of this estimate can be 
reduced to obtaining a suitable estimate on the difference 

|det(/ -Bo + Pae'^PbC-'^Bo) - det(/ - Bo + PaBo)\ 

for b > a > —M. An important technical problem is that the kernels appearing inside 
these determinants are not trace class, so one needs to conjugate appropriately. We refer 
to [QR12a] for the details. The argument for Airy2 is similar. 

As we mentioned, the Airy processes are expected to look locally like a Brownian motion. 
In this direction, it can be shown using the boundary value kernel formulas that the finite 
dimensional distributions of the Airyi process converge under diffusive scaling to those 
of a Brownian motion. The same result was proved earlier by Hagg [Hag08] for Airy2 
using different techinques. In fact, for Airy2 a stronger statement is now available (Corwin 
and Hammond [CHll]), namely that it is locally absolutely continuous with respect to 
Brownian motion. 

Theorem 4.5 ([QR12a], Theorem 3). For any fixed s G M, let B^{-) be defined by Bs{t) = 
£^'^^'^{Ai{s + et) — Ai{s)), t > 0. Then i?e(-) converges to Brownian motion in the sense 
of convergence of finite dimensional distributions. The same holds for B^{-) defined by 

Be{t) = Be{-t), t > 0. 

The proof of this result follows from an explicit computation of 

¥{Ai{eti) <x + y/eyi,...,Ai{etn) < x + yJeyn\Ai{Q) = x) 

and its limit as e — t- 0, see [QR12a] for the details. The same proof works for the Airy2 
process and, in view of (3.4), it should be simple to adapt it to the Airy2-s.i process. 
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4.4. Marginals of Airy2_i.i. The last application of the results of Section 3.2 that we 
will discuss is a proof of the conjecture (1.34) that the marginals of the Airy2->i process 
can be obtained from a variational problem for .42 (t) — on a half-line. The result is the 
following: 

Theorem 4.6 ([QR12b], Theorem 1). Fix a G M. For every m G M, 

sup {A2{t) -t^) <m- min{0, a}^ ) = P(^2^i(a) < m) . 
\t<a J 

The right hand side can be expressed in terms of a Fredholm determinant. Define the 
crossover distributions G"^^^ , for a G M, as 

G2-^i(m) =P(^2^i(a) <m). 

We claim that 

(4.17) G^^\m) = det(l - PmK^Pm), 

where = + and the kernels and K'^ are given by 

Kl{x,y)= I (iAe^"^Ai(x- A + max{0,a}^)Ai(2/ + A + max{0,a}^) 
Jo 



and 



Kl{x,y)= I (iAAi(x + A + max{0,a}^)Ai(y + A + max{0,a}^). 



As noted in Appendix A of [BFSOSb], the kernel i^f^i defined in (1.24) can be expressed 
in terms of Airy functions^: 



KT^^{s, t; X, y) = Lois, x; t, y) + e^tV^-2sVWy-s,^L^ ^ ^^^(^^ ^. 



where 



Lo{s,x;t,y) = -e(-*)^(x,y) = - ^ e-^^-v^ 

POO 

Li{s,x;t,y) = / dA e^("+*) Ai(x - A) Ai(y + A), 
Jo 

POO 

L2{s,x,t,y)= / dAe^(*-'^) Ai(x + A)Ai(y + A) 
Jo 

with X = X — s^ls<o, y = y — t^lt<o, X = X + s^ls>o and y = y + t^lt>o- Using this for 
s = t = a it is straightforward to check that i^'|^|(t, •; t, •) is just a conjugation of the 
kernel K^, and (4.17) follows. 

The fact that crosses over between the GUE and GOE distributions is of course a 

particular case of the crossover property of the Airy2_5.i process, but can be easily obtained 
from (4.17) as well (see the discussion after Theorem 1 in [QR12b]). 

The proof of Theorem 4.6 is similar to (and, in fact, somewhat simpler than) the proof 
of (1.25). Basically, one applies Theorem 3.1 and the cyclic property of determinants to 
compute the desired probability as 

lim ¥{A2{t) < g{t) for t € [-L, a]) = lim detfl - Km + e^^+^^^i^AiOfr „ii^Ai 



L— >c 

with g{t) = t^ + rn and m = m — min{0,Q;}^. An argument similar to the one used in 
Section 4.1 (applying the Baker-Campbell-Hausdorff formula and later checking the result 



^This corresponds to a minor correction of the formula appearing in [BFSOSb], where the exponential 
prefactor appears in front of Lq instead of ii + L2. 
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rigorously, plus some asymptotic analysis to show that an error term goes to in trace 
class norm as L — )• oo) yields 

(4.18) P(sup(^2(t) -i^) <m-mm{0,a}A 
\t<a J 

= det^I - KAiPfn+a^KAi - i^Aie°'^£'m+a2e""^P5R+Q,2i^Ai) • 

Since Kxi = BqPqBq and Bq = I we have by the cyclic property of determinants that the 
right hand side of (4.18) equals 

det(l - PoBoPrn+a^BoPo - PoBoc''^ Qrn+a^e-''^ Pfn+a^ BqPo) • 

Shifting the variables in the last determinant by — m we deduce that 



where 



and 



sup {A2{t) -t^)<m) = det(/ - PmEiP^ - PmE2Pm) , 

t>a 



dX Ai(x - m + 2m + - A)e-2(^-™-"')« Ai(y - m + A) 

-oo 



poo 

E2{x,y) = / dA Ai(x — m + A) Ai(y — m + A). 

Shifting A by m+a^ in both integrals and changing A to —A shows that Ei{x, y) = K^{y, x) 
and E2 = -f^Q, whence the equality in Theorem 4.6 follows since E^ = and E2 = K'^. 
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